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Abstract 

We consider a processor sharing storage allocation model, which has 
m primary holding spaces and infinitely many secondary ones, and 
a single processor servicing the stored items (customers). All of the 
spaces are numbered and ordered. An arriving customer takes the 
lowest available space. We define the traffic intensity /? to be A//i 
where A is the customers' arrival rate and fi is the service rate of the 
processor. We study the joint probability distribution of the numbers 
of occupied primary and secondary spaces. We study the problem 
in two asymptotic limits: (1) m ^ cxd with a fixed p < 1, and (2) 
yO I 1, m —>■ oo with m(l — p) = 0(1). 



1 Introduction 

We consider the following storage allocation model. Suppose that near a 
restaurant there are m primary parking spaces and across the street there 
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are infinitely many additional ones. However, the restaurant has only one 
waiter who serves all of the customers. All of the parking spaces are num- 
bered and ordered; the one with rank = 1 is closest to the restaurant and 
the primary spaces are numbered {1, 2, 3, , m}. We assume the follow- 
ing: (1) customers arrive according to a Poisson process with rate A, (2) the 
waiter works at rate fi, (3) an arriving car parks in the lowest-numbered 
available space, and (4) if there are N customers in the restaurant, the 
waiter serves each customer at the rate fi/N. This corresponds to a pro- 
cessor sharing (PS) service discipline. 

Dynamic storage allocation and the fragmentation of computer mem- 
ory are among the many applications of this model. We define A^^i to be 
the number of occupied primary spaces and N2 to be the number of occu- 
pied secondary spaces. Then we define S to be the set of the indices of the 
occupied spaces, and the "wasted spaces" W are defined as the difference 
between the largest index of the occupied spaces (Max S) and the total 
number of occupied spaces (|S| = A^^i -I- A^2)- Coffman, Flatto, and Leighton 
II2I showed that for the processor-sharing model 



where E\yV] is the expected value of the wasted spaces. Here E[W] = 
6(/(p)) means that there exist positive constants c, c' such that c'/(p) < 
E[W] < cf{p). Also when p ^ I (the heavy traffic case) Coffman and 
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Mitrani obtained upper and lower bounds on E[W] in the form 



2V 1-p ~ ' - 1-p V 6 

A related model, the M/M/oo queue with ranked servers, has been 
studied by many authors d, 0, H, EOl, (13, (III. This differs from the 



current model in that if there are a total oi N = Ni + N2 spaces occupied, 
the total service rate is ^N, as each customer in the restaurant is served at 
rate p. For this model various asymptotic studies appear in 111, m, ffTTfl . In 
particular, Aldous 111 showed that the mean number of the wasted spaces 



is E^] ~ v2ploglogp as p = X/ fi ^ 00. 

A simple derivation of the exact joint distribution of finding A^i (resp.. 



A^2) occupied primary (resp., secondary) spaces appears in 1131 and de- 
tailed asymptotic results for this joint distribution appear in |[9l, while 
the distribution of Max S is analyzed in [7J. In [8J Knessl showed how 
to obtain asymptotic results for the infinite server model directly from the 
basic difference equation. Since the present processor sharing model does 
not seem amenable to exact solution, we shall employ such a direct asymp- 
totic approach here. 

In this paper, we study the joint probability distribution of the numbers 
of occupied spaces in the PS model, letting 7i{k, r) = Prob[Ni = k, N2 = r] 



in the steady state. In part I ||15l we obtained exact solutions for m = 1 and 
m = 2, and developed a semi-numerical semi-analytic method for general 
m. We also derived asymptotic results in the heavy traffic case p t 1/ but 
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with m — 0{1). Here we shall obtain asymptotic results for m ^ oo, for 
the cases < p < 1 and 1 — p = 0{m~^). 

The paper is organized as follows. In section 2 we state the problem 
and obtain the basic difference equations for ^{k, r). In section 3 we sum- 
marize our main results. In section 4 we consider m ^ oo with a fixed p 
(0 < p < 1) and study 7r{k, r) for various ranges of {k, r). In section 5 we 
consider the double limit m ^ oo and p t 1/ with m(l — p) — 0{1). Some 
numerical studies and comparisons appear in section 6. 

2 Statement of the problem 

We let Ni{t) (resp., N2{t)) denote the number of primary (resp., secondary) 
spaces occupied at time t. The joint steady state distribution function is 

Tr{k,r) = TT{k,r;m) = lim Prob[Ni{t) = k,N2{t) = r], < k < m, r > 0. 

t—*oo 

Let p = X/phe the traffic intensity and we assume the stability condi- 
tion p < 1. The pair (A^i, A^2) forms a Markov chain whose transition rates 
are sketched in Fig. 1. The state space is the lattice strip {{k,r) : < k < 
m, r > 0} and the balance equations are 

A; + 1 

(l-^[fc+r>o] + p) 7r(A;, r) = p 7r(A; - 1, r)/[fe>i] + n{k + 1, r)/[fe<„] 

/c + r + i 

r + 1 

+ 1—, — r 7r(/c, r + 1) + p 7r(m, r - l)I[k=m, r>i] ■ (2.1) 

k + r + 1 L - J 
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Here / is an indicator function. The normalization condition is 

oo m 

5^5^7r(A:,r) = l. (2.2) 



r=0 fc=0 



From our viewpoint we will need to consider explicitly the boundary 
conditions inherent in (|2.1[) , so we rewrite the main equation as 



k -\- \ r -\- 1 

(1+p) iT{k,r) = pTT{k-l,r) + — — - 7r(/c + l,r) + — — - 7r(fc,r + l), 

/e + r + 1 /c + r + 1 

0<A;<m, r>0, A; + r>0. (2.3) 



and the boundary condition at A; = m is 

r + 1 

(1 + p) 7r(m, r) = p 7r(m — l,r)H 7r(m, r + l)+p 7r(m, r — 1), r > 1. 

m + r + 1 

(2.4) 

There are also the two comer conditions 

p7r(0,0) =7r(l,0) + 7r(0,l) (2.5) 

and 

(1 + p)7r(m, 0) = p7r(m -1,0) + -^—7T{m, 1). (2.6) 

m + 1 



In (|2.3|) when k = we interpret 7r(— 1, r) as 0. The boundary condition at 
k = min (|2.4|) can be replaced by the artificial boundary condition 

— 7i(m + 1, r) = p 7r(m, r — 1). (2.7) 

m + r + 1 

This is obtained by extending (|2.3|) to hold also at k = m and comparing 
this to (IZ4l) . 



We note that the total number A^i + N2 behaves as the number of cus- 
tomers in the M/M/1 — PS queue, which is well known to follow a geo- 
metric distribution. Thus we have 

vr(A;,r) = (l-p)p^ iV > (2.8) 

k+r=N 

and we can rewrite this as 

N 

Y-K{N-r,r) = (l-p)p^, 0<A^<m, (2.9) 

r=0 

N m 

7r{N-r,r) = 7r(fc, - fc) = (1 - p)p^, > m. (2.10) 

r=N—m k=0 

These identities will provide a useful check on the calculations that follow. 
Fig, 1 A sketch of the transition rates. 
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3 Summary 



Since the analysis shall become quite involved and technical, we collect 
here some of the main results. 

We take m ^ oo, jfirst assuming that p < 1. We use the scaled variables 
X = k/m and y = r/m, so that x is the fraction of primary spaces that are 
utilized. The main approximation we obtain is then given by 

7r{k, r) ~ K{x, y)e'^'t>{x,y). o<a;<l, |/>0 (3.1) 

where 



(t){x, y) = xcpcc + y(t>y + log p - log(s + 1) , 
's + 1 - p - ps e^^'P^^' 



(t>x = log 



s + l- p- s e(i-/')* 

p(s + 1 - p - ps e(^-^)*) 
_p{s + 1 - p) - ps e(i-p)* + (1 - p)(s + 1 - p)e* 

and [x, y) are related to (s, t) via the mapping 

X — x[s,t) 



(3.2) 
(3.3) 

, (3.4) 



[s + l-p-s e^^-"^*] 



(. + l)(l-p)2 
X [(s + 1 - p) e-^^-")* - sp2 - P(l - P)se-*] , 
s 



(3.5) 



(s + i)(i-p) 

Then -^^(a;, y) is given by 



[(1 - p)(s + 1 - p) - ps e-f^ + p(s + 1 - p) e-*] . 

(3.6) 



K{x,y) 



(1 - p)(s + 1 - p)3[s + 1 - p - sp e(i-/')*]3 



(s + l-p-s e(i-/')*)[p(s + 1 - p) - ps e(i-/')* + (1 - p)(s + 1 - p)e*] 



X 



(3.7) 



where j — xtVs — XsUt is the Jacobian of the transformation above. This 
approximation to 7r(/c, r) is explicit in terms of {s,t) but implicit in terms 
of {x, y). However, ior x ^ 1 it becomes much more explicit, with 



7r(/c, r) ~ (1 — p) 



y + 1 



(3.8) 



The above holds for y > and k = m — 0{1), which corresponds to all but 
a few occupied primary spaces. 

Different expansions must be constructed in various boundary and cor- 
ner regions of the strip {{x,y) : < x < 1, y > 0}. For k — 0(1) 
{x — 0{m~^)), which means that only a few primary spaces are occupied. 



(3.9) 



where 

0(0, y) 



yiogp + 



i-p 



log 



y+ v/(l/-l + p)2 + %-(!- p) 
y + - 1 + py + % + (1 - p) 



log y + V(y - 1 + p)2 + 4y + 1 + p + log 2 + logp, (3.10) 



So{y) 



logp + 



log 



y + Vd/ - 1 + py + % - (1 - p) 



2^(l-p)(s, + l-p)3/2y5; 



(s, + l)V(s* + l)2-p 

1 
2 



(3.11) 
(3.12) 

(3.13) 



For r = 0(1) and < a; < 1, which corresponds to having a few secondary 
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spaces occupied and a fraction of the primary ones, we obtain 

Hk,r) ~ i- ' r>l, (3.14) 

.fcO)-(l-p)/ ~ - \_7+p.W-P) ■ ("5) 

Here we wrote the result for r = so as to estimate the deviation of 7r(A;, 0) 
from the geometric distribution (1 — = (1 — 

Near the corner {x, y) = (0, 0) we use the original discrete variables 
(A;, r) to find that 



7i{k,r) ~ p' 



,m+r^-r/(l-p)^^ - p)2-'-/(l-P) T (^1 + 



" 2^/ (3.16) 



and 

^{K 0) - (1 - pjp'^ - -p-m-Va-p) 1^ f ^f^) (3.17) 

p{k,i) = p(i-p)^-^/(^-)r(^) 

X TT^ <P 77^ ^^JN^ , dz. (3.18) 

Here the integrals are over a small loop about z = 0, and these contour 
integrals may be expressed in terms of hypergeometric functions. Near 
the other corner {x,y) — (1, 0) we use the variables n — m — k and r, and 
obtain 

7r(fc,r) ~ r>l (3.19) 

7r(A;,0) ^ (3.20) 



We note that for p < I and m ^ oo most of the probability mass occurs in 
the range k = 0(1) and r = 0, and 7r{k, r) is exponentially small in all of 
the other ranges. 

Defining the marginal distribution by 

oo 

A^(A;) = ^7r(A:,r), < k < m (3.21) 

and 

m 

Af(r) = J2 '^{k, r), r > 0, (3.22) 

fc=0 

we can easily obtain their expansions from the results for n(k, r). For the 
distribution of the number of occupied primary spaces we have 

„m+l p/(l-p) 

M{k)-{l-p)p'' ~ ^ , 0<x<l, (3.23) 

(1 — p)m 

M{k) - {1- p)p'' ~ p™m-i/(i"^)[P(fc,0) + P(A;,l)], A; = 0(1), 

(3.24) 

where P{k, 0) is given by (|4.85|) . For the distribution of the number of 
occupied secondary spaces we obtain 

Ar(r) ~ (1 - p)p-fe+i) (^^±1^^ , ^>o, (3.25) 
Af{r) ~ m{l-pfp^^^- du, r>l, (3.26) 

Jo (1 - p + puy^ 

1-Af(0) ~ m(l-p)V+W du. (3.27) 

Jo I- P + pu 

In particular the mean number of occupied secondary spaces is 

f^rAfir) ~ ( mp"''+\ (3.28) 

^ V2-p; 
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Finally we consider the double limit where m ^ oo and p t 1- We 
introduce the parameter a — m{l — p) = 0(1). On the (x, y) scale we find 
that 

7r(A;, r) ~ m-^/C(x, y)e"**(^'f) ;0<x<l, y>0 (3.29) 

where 

* = + - log(s + 1) (3.30) 

with 

'' + ^-^'- (3.31) 

(3.32) 

(3.33) 
(3.34) 









= log(^ 




1 

s + 1 




s 

s + 1 



1 - si 
1 + s - st 
1 — si + se* , 

(l-si)(l + 2s-si-se~*), 



-[s+{l- st)e-'] 



and 

/C = 



f/2 



v/s(s + 1)(1 - st)(se* - si + 1) V2s(s + 2) - [s(s + 2)^2 - 2i + s(s + 2) - Ije" 



X exp 



ast 



- a{s + 1) 



2(s + l) 

For k — m — 0(1) {x — 1 — 0{m~^)) the expression simplifies to 

7r(A;, r) ~ -^^(t/ + if-m^-<y+^)^ y > 
m y + 1 



(3.35) 



(3.36) 



which shows that the total probability mass in the range k — m — 0{l) and 
y > is asymptotically e^". The remaining mass occurs in the range r = 
and < X < 1, as the results below show that Y^'^=q ^(^) 0) ~ 1 — e~". 
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For k = 0(1) {x = 0{m-^)) we find that 



(3.37) 



where 



*(0,y) 
So{y) 



y + Vy' + % + 2 

2 



+ log 



3/2 



2/+ \/pT4y' 

1 /2 

(y + Vy'' + % + 2) (2/ + V?/' + % + 4) 



(3.38) 



(3.39) 



X exp 



-a I h 



1 



. (3.40) 



2 ■ y+ V2/' + 4y + 2^ 

If y > the analysis and results for 1 — p = 0{m^^) are similar to the 
case p <1. However this is not the case if y = o(l). We give below various 
results for r = {y — 0{m~^)) and r = 0{^/rn) {y — 0(m~^/^)). First, 
for r = 0(1) and < a; < 1 we have 

7r{k,r) ~ - ^'-/2-V-r/2^r/2+l/4g-a(.+l)/2g-2V^(l-V5)^ ^ > (341) 

2 



and 



^ ^-3/4 

A; . . ""^ „-a(a;+l)/2 -2V^(1-Vi) 



7r(/c, 0) - (1 - p)p- ~ ■ -'-e 

2 y/m 



(3.42) 



Note also that with this scaling of p 



(l-p)/ = 



ae 



m 



2m 



+ 0{m-^) 



(3.43) 
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When k,r — 0{1) we have, for r > 1, 



-fc-i 



e''^^-'Uz, (3.44) 



and 



7r(A;, 0) - (1 - p)p^ ~ -av^e-("+^)/'m-5/^e-'^ 



^-fcgl/(l-2) 



(3.45) 

The contour integral may be expressed in terms of a confluent hypergeo- 
metric function. 

The analysis near the corner {x^y) = (1, 0) is much more complicated in 
the case 1 — p — 0{m~^) than when p < 1. We have obtained some partial 
results in the corner range, and there are several nested comer layers that 
must be considered. First when 1 — x = 0(m~^/^) and y = 0(m~^/^) we let 
k — m — ^/rn^ and r = ^/rnR. We find that, for ^,R> 0, 



7r{k, r) ~ m-^/^n{C, R) = m-^/^R-^V{C, R) 



(3.46) 



where 'D{^, R) is expressible in terms of 'D{0, R) via the integral 



1 



exp 



X 



+ 



41og(x/i?) 

{x-R)v{o,x) 



X 



Vlog(x/i?) ' 2x[log(x/i?)]3/2 xV^og{x/R) 



+ 



dx 



no,x) 



2^ J J, 2x[log{x/R)f/' 



exp 



ix-R + O'' 

41og(x/i?) 



c^x (3.47) 
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and V{0,R) satisfies the integral equation 



V{0,R) 



+ 



{x-R)v{o,x) v^{o,x) 



^■Ir I xW^og{x/R) 2x[log(x/i?)]3/2 xV^og{x/R) 



xexp - 7 / dx- (3.48) 

L 4:log{x/R)\ 

We establish several asymptotic results for 'D{^,R) in subsection 5.4, and 
also characterize I> as the solution to a heat equation with a moving bound- 
ary. In subsection 5.5 we consider the scale k — m — 0(m^/^) and r — 
0{y/m) and obtain explicit expressions for 7r{k,r), but these follow simply 
by expanding the result on the {x, y) scale as {x,y) (1, 0) along certain 
parabolic paths. In subsection 5.6 we consider x ~ 1 and r = 0(1), and 

find that with k = m — ^(log m) ^^/m and r > 1 
1 



7r(A;, r) 
7r(A;, r) 
7r(/c, r) 



2 



27r myTogm 



(3.49) 

f+^V«^ o<e<V^, (3.50) 



ae 



2V27r 



^-r/2-l^-r/2g-®VrlogTO 



■00 



= (e-v^)v/b^ = 0(l). 



When r = we obtain 



^{k, 0) - (1 - p)p^ ~ 



7r(A;,0) e ~ ^= 

m v27rmvlogm 



P/2 



(3.51) 



(3.52) 



^(n-a^(-ae^^^ o<e<i 



7r(fc,0)--e- 



2V27r 



J — t 



(3.53) 



0(1) 



(3.54) 



14 



Note that x = 1 — ^(log m)m so that e ~ e [l + a^(log m) / . 



In subsection 5.6 we shall also discuss the scales k = m — 0{^/m\ogm), 
k = m — 0{^/rn), and k = m — 0(1), but we have not been able to resolve 
completely all of the corner layer(s) near (x, y) = (1, 0). Furthermore, the 
analysis also suggests that we may get different asymptotics for the limit 
m ^ oo with 1 — p = 0(m^^/^), which we do not consider here. 



4 Asymptotic expansions for m ^ oo with fixed 

0< p<l 



In this section we examine the case when m ^ oo and p is fixed (0 < p < 1). 
We set 

6 = — , X = — = 6k, y = — = 6r (4.1) 

mm m 

and then let 



7i{k, r) = K{x, y; 6) exp 
For < X < 1, and y > 0, (|2.3)) becomes 



6 



(f){x,y) 



{l+p){x + y + 5) K{x, y; 5) exp 



^4>{x,y) 



p{x + y + S) K{x — 6, y; S) exp 



-0(x - 6,y) 



-(j){x + 6,y) 



+ {x + S) K{x + S, y; S) exp 

+ {y + 6) K{x, y + 5;6) exp 



(4.2) 



-(j){x,y + 6) 



(4.3) 
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The boundary condition (|2.7[) at x = 1 {k = m) becomes 



pK{l,y-5;5)exp 



1 + 5 
l + y + 5 



K{l+5, y; 5) exp 



^<P{l + S,y) 
(4.4) 



We assume that K{x, y; 6) has an asymptotic expansion in powers of 6, 
with 

K{x, y; 6) = K{x, y) + 5K'^^\x, y) + 0(5^). (4.5) 

To compute the leading term for Tx{k,, r) we must compute 0(2;, y) and the 
leading term in (|4.5|) . We divide (|4.3|) by i^(a;, exp[0(a;, and let 5 — *• 
0. Then we similarly divide (|4.4[) by K{l,y) exp[0(l, and let 6^0. 
This leads to the following 'eikonal' equation for 0, for < x < 1 and 

y>o, 

(1 + p){x + y) - p{x + ?/)e"'^" - xe'l'^ - ye"^^ = 0, (4.6) 
and at X = 1 we have the boundary condition 



(4.7) 



We solve (|4.6|) and (|4.7|| for ?/) by using the method of characteristics. 



We write (|46l) as ?/ 



5 V; Y^X5 fy ) 



where 



F = (1 + + - p(x + ?/)e"'^" - xe"^ 



ye- 



(4.8) 
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The characteristic equations for this nonlinear PDE are [H 

dF 



X 



y = ^ 



dx 

dt d(j)x 

dy _ dF 

dt dd),, 



p{x + y)e 



xe^ 



dt 



X + (j)yy, 

-{1 + p) + p e-'^- + e'^ 



dF 

dx 
dF 

-^ = -(l + p)+pe 



(4.9) 
(4.10) 

(4.11) 
(4.12) 
(4.13) 



To solve this system we use rays starting from {x,y) = at t = 0. 



Thus the 'initial manifold' is a; = 1 where the boundary condition in (|4.7|) 
applies, and all the rays start from x = 1 at t = 0. We view x and y as 
functions of s and t. We set 0^ = log T(s, t) and, from (|48l) , (|49l) , (|4.10|) and 
dHH), obtain 

X + y 2p 



T 



-(1 + p), 



x + y 
dT 

— = (T-p)(T-l). 



Solving (|4.15|) gives 



T 



P 



(4.14) 
(4.15) 

(4.16) 



and using (|4.16|) we solve (|4.14[) to obtain 



a; + 7/ = Ci(s)[Co(5)e(^-^)*-p]Vi-^)*. 



(4.17) 



The functions Co(s) and Ci{s) will be determined shortly. 
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From (|4.16|) and (|4.17|) , (|49|) can be written as 



dx 



_Co(s)e(/'-i)* - 1 
Solving (|4.18|l we obtain x{s,t) as 



X = pC^is)[Cois)e^''-'^' ~ p][Cois) - e^^"")*]. (4.18) 



's,t) = [Cois)C,is) - p2Ci(s)e(i-^)* + C2(s)e-''*][Co(s)e(^-^)* - 1], (4.19) 



and then y{s, t) follows from (|4.17|) as 

t) = Co{s)Ci{s){l - pf + C2(s)[e-^* - Co(s)e-*]. (4.20) 

Applying the initial conditions {x{s, 0) = 1, y{s, 0) = s) we find that 

s + 1 



C2{S) 



[Co{s)-pr 

1 Co(s)-p2 



Co{s)-l [Co{s)~p]' 



(s + l). 



(4.21) 
(4.22) 



To determine Co{s) we set (p^ = A{s) and 0^ = B{s) att = 0, and from 
(|4.6|) and (|4.7|) obtain the following equations along the initial manifold 



(l+p)(l + s) = p(l + s)e-^ + e^ + se^, (4.23) 
p(l + s)e-^ = e^. (4.24) 



Solving (523]) and leads to 



=log(s + l), B{s)=\ogp. 



Since e''^« = — y/?/ from (|4.10|) , 0^(1, s) = log[— ?/(s, 0)/s] = B{s) = log p. 
Thus we obtain Co(s) = (s + l-p)/s and then, from (I4.19|I - (I4.22|) , obtain 
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(I331l and ^M- Solving (I4l2l) and (I4l3b then yields ([331) and dM]). Finally, 



we obtain ?/) by integrating (|4.11|) with respect to t. Using integration 
by parts, we can write 

= y {(j)xX + (j)yy)dt = (j)^x + (j)yy - j {(p^x + (j)yy)dt. (4.25) 



The integrand in the last integral in (|4.25|) is —F = 0. We can also check 



from (|3.3|l - (|3.6|) that indeed 4)x/4>y = —y/x. Therefore, (t){x,y) = (p^^x + 
4>yy + f{s) where f{s) is some function of s. Ait = 0, 0(1, s) = \og{s + 1) + 
s log p + f{s) and since ?/ = s at t = 0, 

90(1, 1 



9s s + 1 



+ lOgP + f{s) = (f)y{l, S) = lOgP 



Thus f{s) = — log(s + 1) + constant. We will later show, by asymptotic 
matching, that near the comer (x, y) = (1, 0) the solution Tc{k, r) must be 
0(p'") = O (exp[(logp)/5]). Thus 0(1,0) = logp and we then obtain (Q. 



The value of the constant can also be ultimately found by normalization 
or by using (|2.8|) . Note, however, that our expansion (|4.2[) applies only in 
the domain 0<x<l, y > 0. We shall construct different expansions for 
X ^ and ?/ ~ 0, and also near the two comers (x, y) = {0, 0) and (1,0). 

We can plot the rays, which are given by (|3.5|) and (|3.6|) in parametric 
forms. Each fixed value of s corresponds to a particular ray parameterized 
by t. We sketch the rays in Fig. 2 and Fig. 3. In Fig. 2 we plot the rays for 

s > and t in the range < t < tmax = log[l + (1 - p)/s]/(l - p). At t = 
we have x = 1 and at t = tmax{s) the rays reach x = 0. We also note that, 
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in view of (|3.3|) , (p^ develops a singularity ait = tmaxi-s). However, using 
(|3.5|l and (|3.6|) we can continue the rays for t > tmax, and this continuation 
is sketched in Fig. 3. We see that each ray reaches a minimum value of x, 
where it has a cusp. The cusp occurs when x = y = 0, which happens at 
t = = \og[{s + 1 - p)/(ps)]/(l - p) > t^axi-s). From (HD and dSD we 
also see that when t = tc, x + y = 0, so that the line y = —x is the locus 
of the cusp points. For t > we again have x > and the rays eventually 
re-enter the domain x > 0, and then x increases past x = I. The full rays 
are sketched in Fig. 3, up to the time they return to a; = 1. The figure 
clearly shows the cusps and their locus. However, since the cusps occur 
outside of the domain of interest {0 < x < 1, y > 0} they do not affect 
the solution very much, and we only consider the rays for < t < tmax- 
The origin (0, 0) is the only cusp point in the domain, and this will affect 
significantly the asymptotic structure of 7r(A;,r) for the scale k,r = 0(1). 
We analyze this range in subsection 4.4. 
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-2-10 1 



Fig. 3 A sketch of the extended rays from x — 1. 
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We next calculate K{x, y) in (|4.5|) . Expanding (|4.3|) as 5 — and using 
the fact that satisfies (|4.6|) we obtain the following 'transport' equation. 

""^^xx + pe-^^ + + l) e^^ + + l) e'i'y - (1 + p) 

= [p{x + y)e"'^- - xe-^-j /s:^ - ye'^^ JsTj,. (4.26) 

We note that the right hand side is xK^ + yKy = K, which is the directional 
derivative of K along a ray. The factor that multiplies K in the left hand 
side of (|4.26)) may be rewritten as 



Defining t) to be XtVs — XsUt, which is the Jacobian of the transformation 
from (x, y) to (s, t) coordinates, we find that 



and 



dx_ _ dx^ _ _ XttVs - XtsVt 

OX ox J 

dy _ dyt _ _ ytsXt - yttXs 

dy dy ■ j 



dj 

— = J = xttys + ytsXt - xtsyt - yttXs- 
ot 



Thus, from (|4l2)) , (l4l3l) , and (|4!27l) - (l430l) , we can write (14261) as 



--^-^^x + 0j;) + -Pe- 



2j 2^' 2 
Integrating (|4.31|) using (|3.3|) for we obtain 

, ^ , , n(3+i-p)e-(^-^)^-.pn 

K{x,y) = Ko{s) exp 



(4.28) 
(4.29) 

(4.30) 
(4.31) 



')x + (pv+t 



(4.32) 
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where Kq(s) is a function of s. 



To determine Ko{s) we use the boundary condition at x = 1. From (|4.4|) 
we obtained (|4.7|) to leading order in 5, and at the next order we obtain 



'TO 



Since -s) = ^/(p — s — 1) and (pyyil, s) = 0, (|4.33)) becomes 



+ l)[i^.(l, s) + K,(l, s)] + s) = 0. 



Using the relations 



2(p-s-l) 



Kttj, + KsSy = -Kt— + Ks- = K 
J J 



p — s — 1 



j(s,0) = p-l-s, 



(4.33) 



(4.34) 



(|4.34|) becomes 



Kf , . , OS s 

— + p-1 s + 1 — + -^- = 0, t = 0. 

K ^'^ ' K 1 + s 2 



Evaluating (|4.26|) at t = yields 



Kt s{p + s + l) p 

^ H 7 + s 



(4.35) 



(4.36) 



2(p-s-l) s + 1 
and thus, from (|4.35|) and (|4.36|) , we obtain the following ODE for K at 
t = with respect to s. 



1 



s + l- p s + 1 
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(4.37) 



Integrating (|4.37|) gives 



K{l,s)=C*il 



(4.38) 



where C* is a constant. Equating (|4.32|) at t = with (|4.38|) we obtain -ft'ol-s) 



as 



J, + l_p)3/2 1 



To obtain C* we use the fact that A^^i + follows the geometric dis- 
tribution in (|2.8|) . Expanding and K near x = 1, which corresponds to 
t = 0, the approximation ■K{k, r) ~ /^(x, y) exp[0(x, simplifies to 



7r(fc,r) ~ C* 1 



exp — - — log p exp 



X 



S 



iog(z/+i: 



. (4.39) 



Here we used -ft'(a;, ~ K{l,y) and(f){x,y) = (j){l,y) + (j)x{l,y){x — l) + 

If a; + ?/ is constant, we set x + y = z > 1 and use (|2.8|) with N = z/6 = 
0(5^^). We sum over (A;, r) = (m, N — m) , {m — 1, N — m + 1) , (m — 2, N - 
m + 2), ■ ■ which corresponds to (x, ?/) = (1, z — 1), (1 — 5, 2; — 1 + 5), (1 - 
25, 2 — 1 + 2(5), ■ ■ -. Hence asymptotically for 6 = 1/m (|2.8|) becomes 



A:+r=2/(5 



1 - 



p 



2 + 0(5) 

c _ 



^exp (^^+jjlogp exp[-jlog(z + 0(5))] 



j=0 



Therefore, C* = 1 — p and thus we obtain (|3.7[) . We also obtain 7r(A;, r) near 



X = 1 from (|4.39|) as (|3.8|) , which applies for k = m — 0(1) and y > 0. 
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The ray solution 7r{k, r) ~ K{x^ y) exp[m0(x, y)] is valid in the interior 
of the state space and also near x = 1. However, the expansion breaks 
down near the boundaries a; = and ?/ = 0, and near the comer points 
(x, y) = (1, 0) and (0, 0). We proceed to analyze these regions separately 

4.1 Boundary layer near x = 

We observed from (|3.3|) that (px has a logarithmic singularity as x ^ 0. We 
consider the scale k = 0(1) and y > 0, which corresponds to x = 0{5) and 
r = 0(5-1) ^ We set 

"1 



^■0(0,?/) 



Sk{y]^), 



7r{k,r) = 6" ''exp 
and rewrite (|2.3|) as 

(1 + p)[y + 5{k + 1)] e^(0'^)/^5fc(y; 5) = 5p[y + 5{k + l)]e^(o,,)/5^^_^(^. ^) 

+ {k + l) e^(°'^)/^5,+i(y; 5) + (y + 5) e^(°'^+^)/^5,(|/ + 6; 6) (4.40) 

for k > and ?/ > 0. The boundary condition at A; = is 

{l + p){y + 6) e'^^^'y^/'Soiy, 6) = e'^^^'^^/'S^ {y; 6) + {y + <5)e'^(°'^+^)/^^o(z/ + S;6). 

(4.41) 

We expand Sk{y; S) in the form 

Sk{y;S) = S,{y) + 6Si^\y) + 0{6'). 
In (|4.40|) and (|4.41[) we let 5 ^ and obtain the following equation for the 



leading term Sk{y) 



ri + p_e^«(Ml^ = (A:+i)%tii^, k>0. (4.42) 

Sk{y) 
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The general solution to (|4.42|) is 



and thus 



7i{k,r) ~ S-'^e'^i^'Vys So{y)jT [i + p _ e<^«(o,y)]\ (443) 



It remains to determine the constant u and the function So{y). 

We next asymptotically match (|4.43|) to the ray expansion exp[0(x, 



in an intermediate limit where x ^ but k = x/S ^ 00. Using Stirling's 



formula for k\ the expansion in (|4.43|) for k ^ 00 becomes 



7r(A;,r) ~ S'^e^^^'^^^^ Soiy)^=exp[-klogk - klogS] [l + p ~ e^'^^'^^ y 

V^nk 



— exp 



V2 



TCX 



^ |0(O, y) + x + x\og (J^^ + X log [1 + p - e'^^(°'^)] } 

(4.44) 



Here we rewrote the result in terms of x. 

We expand the ray solution as x 0. Along x = we can explicitly 
invert the transformation from (x, y) to ray coordinates. When x = we 
have t = tmaxis) and s + 1 — p — se^^^''^* = 0. Then from (|3.6|l we find that 
s and y are related by y = s(s + 1 — p)/(s + 1). We thus define = s^{y) 
by (|3.13|) . Thus, a ray that starts from (x, y) = (1, s*) hits the y— axis at the 
point {0,y). Then we use (|3.2|) and (|3.4|) and evaluate these expressions 
at t = tmaa and s = s^,, to obtain explicit expressions (|3.10|) and (|3.11|) for 



0(0, y) and 0j,(O, ?/). We note that 0^(0, y) = log p — tmax{s*)- Also, from p.3|) 
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and (|3.5|) we obtain 



— log x+ log 



[s + l- p- ps e^^-p'^*] [{s + l- p) e-(i-^)* - _ ^(i _ p)se-*] 



(s + l)(l-p)2 



(4.45) 



which is exact for all (x, y) and indicates the logarithmic singularity in 
as X 0. By expanding (|4.45|) for s — and t tmax we obtain 

+ 1 - p)(l + p - pe"*'' 



-logx + log ^ — +o(l) 

- logx + logy + log [l + p - e'^''(°'^)] + oil), 



(4.46) 



which integrates to 

0(x,?/) = 0(O,t/)-xlogx + x + xlog?/ + xlog [l + p - e'^''^°'^^] +o(x). (4.47) 



The above precisely agrees with the exponential terms in (|4.44|) . The ex- 
pansions will thus asymptotically match provided that z/ = — 1/2 and 



So{y) = lim v2nx K{x,y) 



By expanding j = Xtys — x^yt and then K in (|3.7[) for s — and t — > t^ax/ 
we obtain after some calculation 

v/2^(l-p)(s, + l-p)3/2^ 



So{y) 



(S, + I)v/(S* + 1)2-P 

l-p)v/2^ [(1 + p)y + (1 - p) v/(y - 1 + p)2 + 4y + (1 - p)^ 



1/2 



1/2 + 2(1 + p)y + + 1 + p) ^(y - 1 + p)2 + 4y + (1 - p)^ 



nl/2 



(4.48) 



This completes the determination of u and Soiy) in (|4.43|) , and verifies the 
matching between the (x, y) and (A;, y) scales. 
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4.2 Boundary layer near ^ = 

We consider small values of y, and return to the original discrete variable 
r = y/5, thus setting 

7r(A;,r) = n^(x;5) r > 1 (4.49) 

and 

7r(A;,0) = (l-p)p'= + no(x;5). (4.50) 

We assume for now that < x < 1, and will treat the corner regions, 
where x ~ or x ~ 1, separately. Also we expect that for large m and 
p < I, the secondary servers will rarely be needed. Thus we expect that 
7i{k, r) will be mostly concentrated along 7r(k, 0), with this function being 
roughly geometric in p. Thus we wrote 7r(A;, 0) in a form so as to estimate its 
deviation from a geometric distribution. We also note that (1 — p)p''I{r=o} 
is an exact solution to (|2.3|) , and also satisfies the corner condition (|2.5|) at 
{k, r) = (0, 0). It fails to satisfy the full problem in (|2.1|) only due to the 
boundary condition aik = m. With (|4.49|) and (|4.50|) , (|2.3|) becomes 



(l + p)n,(x;5) = pIir{x-5;5) + —!^^^—Iir{x + 5]5) 

X + o[r + 1) 

X + d{r + 1) 

(l + p)no(x;5) = pUo{x-6;6) + Uo{x + 5;5) + -^n,{x;5), r = 0. 

X + d 

(4.52) 

We expand Ilr{x; 5) as 

n,(x; 6) = n,(x) + (5n(^)(x) + 0{6^). 
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From (|4.51|) and (|4.52|) we obtain to leading order the following equations 



x{l - p)n; - rUr + (r + 1)U. 



r+1 



x(i-p)n'o + ni = 0. 



(4.53) 
(4.54) 



Equations (|4.53)) and (|4.54|) express the 'current' value Ilr{x) in terms of the 
'future' value Ilr+i{x). However, for sufficiently large r we can use the ray 
expansion to compute TT{k, r) asymptotically. Expanding K{x, y) exp[m0(x, y)] 
for y ^ and < x < 1 leads to 

7r{k,r) = K{x,0)exp{m[(f){x,0) +y(f)y{x,0) + 0{y^)]} 



Ll p''»|-g</'H{a;,0)jr 



1 - p + pe 
(1 -p)2p"'+^(e"*)'' 



(4.55) 



[1 - p + pe-^Y+^ 

Here we used the fact that y = implies that s = 0, and thus (cf. (|3.5|) ') 
x(0, t) = e^P~^^^ and e^* = x^^^^^''\ The ray that starts from the corner (1, 0) 
is the line segment y = 0, 0<x<l and this ray reaches (x, 0) when 
t = — (logx)/(l — p). Therefore, from (|4.55|) , as ?/ — -K{k,r) becomes 
(|3.14|) . We can easily check that (|3.14|) satisfies (|4.53)) . Thus we have the 
leading term n^(x) for r > 1. To obtain no(x) we solve (|4.54|) with Ili{x) 
computed from the right hand side of (|3.14|) with r = 1. We thus obtain 



no(x) 



p 



c 



1 - p + px^/^^"p\ 

where C* is a constant. Since vr(0,0) = 1 — p, we expect that no(0) = 
p"^[C + {1 — p)] = 0, which implies that C = — (1 — p). Therefore, for 
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r 



we have (|3.15|| . In subsection 4.4 we will consider n{k,r) for r and 



k = 0(1), and will give a more precise argument, based on asymptotic 
matching, to show that no(0) = 0. 

We know that for k < m, J2n=o ^(^ — n,n) = (1 — p)p^ exactly, and this 
implies that Y1^j=q^j{'^ ~ ^J'^ ^) = 0- Since < s < 1, asymptotically we 
should have X]j°=o nj(a;) = 0, which is indeed true since 

^ fl - o)2^m+r r/(l-p) Q _ N m+1 l/(l-p) 

4.3 Corner layer near {x,y) = (1, 0) 

We consider (x, y) near the comer (1, 0). We again use the discrete variable 
r and consider k = m — 0(1) (or 1 — x = 0(5)). Thus we define n and L by 

n = m — k, ii{k, r) = L[n, r; m). 

The main balance equation (|2.3|) , written in terms of (n, r), becomes 

m — ra + l 

(1 + p)L[n, r; m) = pL[n + 1, r; m) H L(r;, — 1, r; m) 

m — n + 1 + r 

r + 1 

H L{n,r + l;m). (4.56) 

m — n + 1 + r 

Here we used 7r(/c±l,r) = L(n=i=l, r; m). The artificial boundary condition 
(|2.7|) becomes 

TTl ~h 1 

L(-l,r;m) = p L(0,r - l;m), r > 2, (4.57) 



m + r + 1 
and the comer condition in (|2.6|) is 



(1 + p)L(0, 0; m) = pL(l, 0; m) + ^^L(0, 1; m). (4.58) 

m + 1 
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By requiring (|4.56|) to hold also at n = 0, thus defining L(— 1, r; m), we can 
use (|4.57[) (or (|2.7[) ) instead of (|2.3|) . Then (|4.58|) may be replaced by 



L(-l,0:m) = 0. 



(4.59) 



The condition at the other corner (0, 0) corresponds asymptotically to n 
oo and will play no role. 

We expand L for m — > oo (or 5^0) with 

L(n, r; m) = L{n, r) H L^^\n, r) + 0(m~^). 

m 

To leading order we obtain from (|4.56|) 



1^1 + p)L{n, r) = pL{n + 1, r) + L(?t, — 1, r), r > 



and (l457b and (l439b lead to 



L(-l,r) = pL(0,r-l), r > 2, 
L(-1,0) = 0. 



(4.60) 

(4.61) 
(4.62) 



We can infer the solution of (|4.60|) - (|4.62[) be expanding the ray solution 



as (x, (1,0). This can be obtained by simply letting x — 1 in (|3.14[) , 
which suggests that 

L(n, r) = (1 - r>l, (4.63) 



which is independent of n. We can easily verify that (|4.63|) satisfies (|4.60|) 
and (|4.61|) . To obtain L{n, 0) we let x ^ 1 in (|3.15|) and recall that k = m — n. 
Hence, 

L(n,0) = (l-p)(p""-p)p'". (4.64) 
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We also note that (|4.58|) implies asymptotically that 



;i + p)L(0,0)=pL(l,0), 



which is indeed satisfied by (|4.64[) . Thus we have obtained the leading 



term for 7r(fc, r) in the comer region in a very simple form, with (|3.19|) and 
(|3.20|) . We note that when r = and k = m — 0(1) even the leading 
term indicates a deviation from the geometric distribution [1 — p)p^ = 
(1 -p)p'"-". 

As a check we verify that (|2.8|) is satisfied asymptotically. In terms of n 
this identity becomes 



l-p)p^= 7r(A;,r)= L{n,r;m). 



(4.65) 



k+r=N m—n+r=N 

We consider N = m + N with m, N ^ oo but with iV = 0(1). Then the 
comer range approximation can be used to approximate 7r(A;, r) in the sum 
in (I4.65D . We may have iV<OoriV>0. IfiV<0, k + r = m + N < m and 
the approximations in (|3.19)) and (|3.20[) yield 



L{n,r;m) = L(r — iV, r; m) ^ L(r — iV, 



r=n+Af 



r=0 



r=0 



L(-iV,0) + J]]L(r-iV,r) 



;i-p)(p^-p) + E(i-/^) 



r=l 
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so that (|2.8|) holds for N < 0. For N > we have N > m and then neces- 
sarily r > 1 so that 

m+N 



L{n,r;m) = ''^^ L{r — N,r;m) 



r=n+N r=N 

oo oo 

~ J]L(r-iV,r)=p-(l-p)2^p'- 

y=JV r=7V 

which again evaluates to (1 — p)p"^~^^ and m + N = N. When N = m 
(N = 0) either of the above calculations apply, since n{m, 0) ~ L(0, 0) can 
be computed by either of the formulas in (|3.19|) and (|3.20|) . This means that 



7r(fc, 0) can be computed by setting r = in the expression for r > 1, but 
only if = m. This observation was also used in the analysis in section 5. 

4.4 Corner layer near (x, y) = (0, 0) 

We consider near (x, y) = (0, 0) and go back to the original discrete vari- 
able {k, r), with k,r = 0(1). We also set 

Tr{k, r) = p^m-'''^^~P~^P{k, r; m), r > 0, (4.66) 

and 

7r(A;, 0) = (1 -p)p'= + p'"m-^/(^-^)p(A;,0;m), r = 0. (4.67) 

We also require that this comer expansion match to the expansion in sub- 
section 4.2, in an intermediate limit where k ^ oo but x = k/m 0. 
By expanding (|3.14|| and (|3.15|l as x — * 0, we see that x'^^^^~p^ becomes 



0(m-''/(^-^)) on the fc-scale. We thus scaled 7r{k,r) to be 0(p'"m-''/(i-'')) 
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expecting that P{k, r; m) will be 0(1). Note also that 7i{k, 1) and 7i{k, 0) — 
{l—p)p^ are scaled to be of the same order in m. Then the balance equation 
(|2.3|) becomes 

k -\- 1 

(1 + p)P(k,r]m) = pP(k — l,r]m) + P(k + l,r]m) 

k + r + 1 

m-^/'^^-p''P{k,r + l]m), k > 0, r > 0, 



k + r + 1 

(4.68) 

[1 + p)P{k,0;m) = pP{k-l,0;m) + P{k + l,0;m) 

+ P{k,l;m), k>0,r = 0. (4.69) 



Expanding P{k, r; m) as P{k, r; m) = P{k, r) + o(l) we obtain from (|4.68[) 
and (|4.69|) the following equations for the leading term 

{l + p)P{k,r) = pP{k-l,r)+ , Pik + l,r), r > 0, (4.70) 

k + r + 1 

{l + p)P{k,0) = pP{k - 1,0) + P{k + 1,0) + P{k,l), k>0,r = 0. 

k + 1 

(4.71) 

Here we define P(— 1, r) = so that (|4.70|) holds for k > 0. However, we 



must now consider the corner condition (|2.5|) . Since 7r(0, 0) = 1 — p, (|4.67|) 
implies that P(0, 0) = and thus (|2.5|) asymptotically becomes 

= P(1,0) + P(0,1). (4.72) 



This is also consistent with (|4.71|) at A; = 0. We shall first analyze (|4.70|) for 



r 



> 0, and then solve (|4!7T|) for P{k, 0). 
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For r > we use the generating function ^{z) = ^ {z; r) = -^(^' ''")^^- 

Multiplying (|4.70|) by z'^ and summing over k > gives 



[-pz' + (1 + p)z - 1] e'(^) = [p(r + 2)z - (1 + p)(r + 1)] ^z)- (4.73) 



The general solution to (|4.73|) is 

e(z;r) = e.(r)(l-;.)-i-'-/(i-^) 
where is a function of r. In particular. 



^ ^ l-pr/(l~p) 



1 — pz 



e(0;r) = Ur)p'~'"'^^'''^ = P(0,r). 
Inverting the generating function, we can write for r > 0, 

where the integral is a complex contour integral along a small loop around 
z = 0. We obtain by matching. Letting x ^ {x = k/m) in (|3.14|) , and 
invoking the asymptotic matching condition shows that 

P{k,r) ~ (1 - p) (p^) k'/^^-p\ k^oo. (4.75) 

We now expand (|4.74|) as /c ^ oo and verify that (|4.75|) is satisfied. This will 



also determine the function ^*{r). As A; ^ oo the behavior of the integral 



in (|4.74|) is determined by the singularity closest to 2; = 0, which occurs at 
z = 1. Thus, expanding the integrand in (|4.74|) near z = 1 we obtain 

l-pr/(l-p) Y(k+1 + YT^) 



P 



^-pJ r (i + j^^) km 



Ur) { 7 ^. (4.76) 

^-PJ kB(k,l + j^^ 
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Expanding the Beta function as — > oo yields 



Bik,i + 1 ~ r ( 1 + ] k-^-^i^^-p\ 



i-pj V i-p 

and thus (|4.76|) simplifies to 

p \ l-pr/(l-p) ^r/(l-p) 



Therefore, (|4.75|) is satisfied if 



p \^-PJ \ 1-p 

We show that (|4.74|) , when expanded for r — > oo asymptotically matches 
to the expansion in subsection 4.1 as y — 0. Letting y — > in (|3.10|) , (|3.11|) , 
and (|4.48|) we obtain 



S,{y) ~ (l-p)3/2v^, 

±(n \ ^ , iM , l/logy 2?/log(l-p) y , 2n 

</' 0,?/ = ?/ + l logp + - + 0y , 

1 -p 1 -p 1 -p 

1 2 



Thus as y — »• (i/ = r /m) the asymptotic behavior of (|4.43|) is 



7r(fc,r) ~ v^(l_p)3/2-2r/(i~p)(i + ^)y+m!^ ry''^^ (4.77) 

This is an approximation to n(k,r) that applies in the matching region 
where k = 0(1), r ^ oo, r/m 0. 



Next we expand the corner approximation (|4.66|) , which is given by 



(|3.16|) , as r — > oo. We expand the integrand in (|3.16|) around ^ = by 
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setting z = u/r and evaluate the integral, using 



2-Ki 



u-^-\^^^P'^''du = —{I + pf . 

K ! 



Using Stirling's formula we can expand the Gamma function in (|3.16|) as 

r/(l-p) 



r 1 + 



27rr 



-r/(l-p) 



OO. 



(4.78) 



1 — p J \ 1 — p \1 — p^ 

Therefore, (|3.16|) as r ^ oo agrees with (|4.77[) . 

Next we verify that the comer approximation (|3.16|) matches to the ray 
solution K{x,y) exp[m(f)(x,y)], in an intermediate limit where A;, r — > cxo 
but X = k/m, y = r/m ^ 0. We expand (|3.16|) for k and r simultaneously 
large, writing the integrand as 

exp[H{z)] _ 1 
z{l-z){l- pz) ^ z{l - z)(\- pz) 

X exp 



pr T 

-k log z + log(l - pz) - logfl - z) 

1- p 1 -p 



and using the saddle point method. The saddle point(s) satisfy H'{z) = 
so that 



k 
r 



+ 



1 



1 — p \ 1 — pz \ — z J 
(|4.79|) defines the saddle as a function oik/r and the saddle zq = ZQ^k/r) is 
given explicitly by 



2p 



1 + P+ Wl -2p 



k — r 
k + r 



(4.80) 



We note that zq as k/r and zq 1 as k/r oo. Evaluating 
the integral along the steepest descent path through the saddle point we 
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obtain 



1 



v/27rF"(^o) {1 - zo){l - pzo) z, 

1 



fc+i 





I- Zq 



(^o) Zo{l - Zo){l - pzo) 



(4.81) 



em {1 — p)^(l — Zq) 

In the ray solution we let x — > and ?/ — > 0, which corresponds to s — > 

and t oo, with a;/?/ fixed (corresponding to se^^^^^* = 0(1)). We find 

that in this limit 

1 - p - se(i-^)* 
1 - p - pse(i-p)* ^ 



logp 



1-p 



[log?/ - 21og(l - p) + plog(l - p2;o) - log(l 



log(l + s) 



~ s ~ 



1-p 



Therefore, as (x, y) — > (0, 0), in terms of zq = zo{x/y) = zo{k/r), 



exp[m(f){x,y)] ~ p'^+'z^'' 



1 r (1 - pzoY 
e m (1 — p)^(l — Zq) 



r/(l-p) 



This agrees precisely with the exponentially varying terms in (|4.81|) . From 
(14.791) and (14.801) we also obtain 

1 + r 



(1 - Zo){l - pZo) 



(1 +p) (k + r) 



2p 
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where 



It follows that 



{l+py{k + r) 
(1 - 



4p X 
{l + p)Hx + yy 



(4.82) 



(4.83) 



F"{zo) Zo{l-Zo){l-pzo) + VV-V2- 

By asymptotic matching, the above should agree with the expansion of 
K{x, y) as x,y ^ 0. As s ^ and t ^ oo with se*^^^'^^* = 0(1), from (|3.7[) 
we obtain 

(1 -p)[l -p-pse(i-^)*] 



K{x,y) 



(4.84) 



[1 - p - se(i-^)*] V2[i - p + pse(i-^)*]V2 • 
Since x(s, t) ~ [1 - p - se^^"^)*] [(1 - p)e-'^^-P^^ - sp^]/{l - pf and ?/(s, t) ~ 
s(l — p) in this limit, (|4.84|) agrees with (|4.83|) . This completes the matching 
verifications. 

It remains to compute P(/c, 0), by solving (|4.71|) and (|4.72[) . We use the 
generating function 

oo 

eo(^) = 5^P(A:,0)^^ 

multiply (|4.71|) by 2*^, and sum over A; > 0, to obtain 

1 



(l+p)z-pz^-l^A; + 



k + l ^ ' ' 



and thus 



P{k,0) = z-^-^Uz)dz 



1 



27ri / (1 --2)(1 - pz) 
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.i=o ^ 



(4.85) 



Evaluating the contour integral leads to 

k-l 



1- pf'-^ p(j, 1) 

j=0 ^ 



p J + 1 



(4.86) 



Finally we verify the matching between (|4.85|) and (|3.15|) . To evaluate 
the integral in (|4.85|) as A; ^ oo, we expand the integrand around z = and 
obtain 



Pik,0) 



1-p 



P{k-l,l) , P(fc-2,1) ^ P(fc-3,1) ^ ^ ^^^^^^ 



A; 



fc- 1 



A;-2 



p ^ (n - l)i/(i-p) 



- p 
p 



i-pji 



n=l 



n 



p) 



(4.87) 



Therefore, from (|4.67[) and (|4.87|| we obtain 



7r(fc,0)~(l-p)p'^-p-+M^ 



i/(i-p) 



(4.88) 



This applies for k oo but with A;/m ^ 0. If we let x — in (|3.15|) , we 
obtain precisely (|4.88|l , which verifies the matching. Thus, for r = the 
approximation to 7r(/c, r) must be computed by using (|4.74|) with r = 1 in 
diiHSl) and then dHZl with P{k, 0; m) ~ P(A;, 0). 
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4.5 Marginal distributions 



We define the marginal distributions by (|3.21|) and (|3.22|) . First, we discuss 



Ai{k) for various ranges of k. Summing (|2.1|) over r, we obtain 

oo oo 

W'>o] +p){k + r + l) 7r{k, r) = pl[k>i] Y,(^ + ^ + ' 1' 

r=0 r=Q 

oo oo 

+ {k + l)J2 ^(^ + 1, ^) + J^'^r + l)7T{k, r + 1). (4.89) 

r=0 r=0 

If we define Mi{k) = J2T=i ^ ^(^' then (|4.89|) can be written as 



+ - A<(/e - 1)] - {k + l)[M{k + 1) - M{k)] 

= 7r(0, 0)/[fc=o] - p[Mi{k) -Mi{k- 1)] (4.90) 

where A^i(-l) = and Mi-1) = 0. Setting Mik) = (1 - p)p^ + M(A;), 
from (|4.90)) , we obtain 



[M(A;+l)-M(fc)]-p[M(A;)-M(A;-l)] = ^ j A^i(fc)-A^i(fc-l)] (4.91) 

(A; + 1) 

where M(-l) = 0. 

On the X scale (corresponding to A; ^ oo, < x < 1), if we set M{k) = 
p'^Mix; 6) and Mi{k) = p™Mi(a;; 6), dHH) becomes 



(x + 6)[M{x + 6;6)~ M{x; 6)] - p{x + 6)[M{x; 6) - M{x - 6; 6)] 

= Sp[Mi {x; 6) -Mi{x-6;6)]. (4.92) 

Expanding M{x; 5) and Mi(x; 5) as 

M{x; 6) = 6M{x) + 0(5^), Mi(x; 5) = Mi{x) + 0(5), 
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from (|4.92|| , we obtain to leading order 

M'ix) = ^MM. (4.93) 

1 — p X 

For < X < 1 7r(fc, r) is maximal in the range r = 0(1) for r > 1. Thus we 
use (|3.14[) to evaluate the sum in (|3.21|) , to obtain 



A<i(A;) = p'"Mi(x;5) 



rn 

r\-/ 



^ [1 - p + pxva-p)] ^ \[i_p + p^i/(i-p)] 



This gives the mean number of occupied secondary spaces, if the number 
of occupied primary spaces is A; = mx. 

To leading order we have M.{k) ~ (1 — p)p^ ~ Tt{k, 0). To obtain the 
correction (second) term for M{k), we need to solve (|4.93|) and then use 



(|4.94)l , which yields 



M(a;) = -i— / -M((?)c/^ = xPl^^~P\ (4.95) 

^ ~ P Jo 1~P 



Therefore, we obtain p.23|) . Using (|3.23|) we note that 



„m+l /•! 

J]A<(fc) ~ 1 -p"^+i + / x^/^^-^^-Mx = 1. (4.96) 

Thus, the marginal is properly normalized up to order 0{p'^). Expression 
(|3.23|) breaks down as x ^ 0, since then we must use the corner expansion 
for 7r(/c, r) to evaluate (|3.21|) . 
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For the k scale (x = 5k, k = 0(1)) we must use the approximation to 
n{k, r) valid for k,r = 0(1). Thus from (|4.66|) and (|4.67[) we obtain 



^m^^i/(i-p) [p^^^ 0) + P(A;,1)], 



M{k) 

Mi{k) ~ p'"m-^/(^-'')p(fc,l) 



(4.97) 
(4.98) 



and we can easily check that (l497t and (l498b satisfy dUT). Then P(A;, 1) 
and P{k, 0) are given by (|4.74)l and (I4.85|) . 

We verify the matching between the k and x scales. For k — > oo, we use 
(|4.98|) and write (|4.91|) asymptotically as 



(1 - p)M'(A;) ~ p-+im-i/a-p)^___^[P(A;, 1) - P(fc - 1, 1)]. (4.99) 



Also for k oo, we use (|4.75|) in (|4.99|) and obtain 

„m+2 



and thus 



M\k) 



M{k) 



P 



(1-p) 



P' 



m+l 



^-l/(l-p)^l/(l-p)-l_ 



(4.100) 



(4.101) 



(1-p) 

This matches to M(k) on the x scale, in view of (|3.23[) . 

Near the boundary of x = 1, to obtain (A;) we use the corner (x, y) = 
(1, 0) solution in subsection 4.3. From (l3l9l) (or irab ) and (|3^ (or (1461) ), 
we obtain 

oo 

Mik) ~ J](l - p)V'"+'- + (1 - P)(P~" - P)P"^ = (1 - P)P'"-" = (1 - P)p'. 

r=l 
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Thus the leading term for the marginal Ai{k) is the geometric distribution 
{1 — p)p^, for all ranges of < /c < m. The correction term is roughly 0{p'^) 
and contains the effects of the (rarely used) secondary servers. 

Next we discuss M{r), the marginal distribution of the number of oc- 
cupied secondary spaces, for various ranges ofy = 6r. On the y scale (r = 
0(m)), we use the fact that 7r(k, r) is maximal in the range k = m — 0(1). 
By summing (|3.8)) over k < m we obtain the distribution of N2 as 

m / _|_ 1 _ \ 

Af{r) = J] 7r(m - n, r) ~ (1 - p) ^ / Z^'" + 1)"^ 

which can be written as (|3.25|) . (|3.25|) applies only for y > 0, and has a 



singularity asy 0. 

On the r scale (r = 0(1)), we integrate (|3.14|) over < x < 1 and obtain 
(13261) from 



Af{r) ~ m(l - p)2p'"+" / :rjr-^^dx. (4.102) 

Here we approximated the sum in (|3.22|) by an integral and changed vari- 
ables using u = x^/*^^ ''^ . To check the matching between the r scale and the 
y scale we let r ^ oo in (|3.26[) . We expand the integrand in (|3.26|) around 
u = 1 since it is an increasing function as r ^ oo. If we set u = 1 — w/r, 
(|3.26)) becomes 

(1 — w /rY~^ dw 



Af{r) ~ m(l - p) V 



3 jm+r 



(1 — pw/rY~^^ r 

1 - p)V"^+^ - e-^'-P^^dw = (1 - p) -.(4.103) 
r Jo r 
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We can easily check that letting y in (|3.25|l and using y = r /m, we also 
obtain (I4l03l) . 

We need a different expression for A/'(r) when r = 0. From (|3.15|) , again 
summing over < A; < m ( and integrating the correction term over < 
X < 1) we obtain 

m 

Ar(0) = 5^7r(A;,0) 



fc=0 

m 



lit It 

k=o -^0 



X 



V{i-p) 



1 — m(l — p) p 



2 



1 ^/l-^ 



1- p + pu 

The normalization condition is satisfied since 



1 - p + pxi/(i-p) 
du + Oip""). 



dx 



(4.104) 



r=l ^ 

= m(l-p)2p"+^ 



u 



- p + pu 



pu 



p + pu 



r=l 

du ~ 1 -7V(0). 



du 



I - p + pu 

We have thus estimated the difference 1 — A/^(0), which gives the exponen- 
tially small probability of having to use the secondary spaces. Note also 
that the mean number of occupied secondary spaces is, from (|3.26|) , 

r{puy 



oo „i oo 

= m(l-p)p™+i / u^-Pdu 



- p + pu)' 



l-p 



u 



2-p 



1 - p + pu) 

m+1 



du 



mp 
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5 Asymptotic expansions for m ^ oo with p t 1 



In this section, we study the heavy traffic case, in which p ] I, with m 
oo. We introduce the parameter a, with 



p = I = 1 — aS, a = 0(1). 



(5.1) 



We shall again analyze (|2.3|) for different ranges of {k,r) (or {x,y)). For 
some of the ranges the analysis will closely parallel that of section 4, which 
had p < 1. However, when y ^ the analysis is much different, and 
we will show that the scalings r = O(v^) {y = 0{\/6)) and r = 0(1) 
(y = 0{5)) will lead to very different types of asymptotics. 

We first consider < x < 1 and y > 0, and set 

'1 



7r{k, r) = V{x, y) = /C(x, y; 6) exp 
Then (|2.3|) becomes 



:^{x,y) 



(2 — a6){x + y + 6) /C(x, y; 6) exp 



-<i/{x,y) 



[1 — a6) {x + y + 6) IC{x — 6,y; 6) exp 



-<i/{x-6,y) 



-^{x + 6,y) 



+ {x + S) JC{x + S, y; 5) exp 

+ {y + 5) /C(x, y + 5;5) exp 



-<i/{x,y + 6) 



(5.2) 



The boundary condition (|2.7[) at x = 1 (A; = m) becomes 



^1— a(5)/C(l, y—S; 6) exp 



-vl>(l,y-5) 



l+y + S 



/C(l+(5, (5) exp 



^^(1 + 5,?/) 
(5.3) 
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We assume that /C(x, S) has an asymptotic expansion in powers of S, with 

JC{x, y; S) = S}C{x, y) + 5^1C'^^\x, y) + O^S"^). (5.4) 

To compute \1/ (a;, y) we divide (|5.2[) by 5K{x,y) exp[\l/(a;, and let5 0. 
Then we similarly divide (|5.3)) by 5/C(l, exp[\E'(l, and let 5 ^ 0. We 
thus obtain the following equation for \E', for < x < 1 and ?/ > 0, 

(x + y)(2-e"*^) -xe*^ -ye*'' =0, (5.5) 

and at x = 1 the boundary condition 

{l + y)e-^y^^^y) =e^^^^^y). (5.6) 



We solve (|5.5|) and (|5.6|) for ^(x, y) by using the method of characteristics. 



as we did in section 4 again writing (|5.5|) as \E', \E'a;, ^y) = where 

J- = (x + (2 - e"*^ ) - xe*^ - ye*^ . (5.7) 

The characteristic equations for this nonlinear PDE are HI 

dx dJ-" 



X 



dt 



X + i/)e-*^ - xe*^ (5.8) 



X 



d-^ 

= -9^ = -2 + e-'^- + e*% (5.11) 
ox 

Fit 

= -^ = -2 + 6-^^^ + 6*^ (5.12) 
dy 
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To solve this system we again use rays starting from (x, y) = (1, s) at t = 0. 
We set = log t(s, t) and, from ^7} - dSH) and (ISlTI) , obtain 



x + y 2(1-1 



T 



The solutions to (|5.14|) and (|5.13|l are 



T 



1 



1 



t + CoisY 
x + y = C,{s)[t + Co{s)-lf. 



Using (|5.15|) and (|5.16|) in (|5^ we obtain 

1 



dx 



t + Co(s) 



x = Ci{s)[t + Co{s)][t + Co{s)-l]. 



The general solution to (|5.17|) is 



= [t + Co(s)]{Ci(s)[t + Co(s) - 2] + C2(s)e-*} 



Subtracting (|5.18|) from (|5.16|) we obtain 



(5.13) 
(5.14) 



(5.15) 
(5.16) 



(5.17) 



(5.18) 



yis,t) = Ciis) - C2{s)[t + Co(s)]e-*. (5.19) 

Applying the initial conditions {x{s, 0) = 1, y{s, 0) = s) we find that 

s + 1 



C2is) 



1 Co(s) - 2 



Co(s) [Co(s) - 1]' 



(. + 1). 



(5.20) 
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To determine Co{s) we set = and = f3{s) at t = 0, and from 
(1531) and dSH) obtain 



2(1 + s) = p(l + s)e-" + + se'', 



(5.21) 
(5.22) 



Solving (15.211) and (|5^ leads to 



a{s) =log(s + 1), Pis) =0. 

Therefore, = a(s) = log(s + 1) = logt = log[l - l/Co(s)] at t = 0. This 
implies that Cq{s) = — 1/s, which then gives Ci{s) = 6*2(5) = s^/(s + 1). 
Thus we obtain x{s, t) and y{s, t) as (|3.33|) and (|3.34|) . From (|5.15|) and (|5.9[) 
we also obtain (|3.31|) and (|3.32|) . In the same manner as in section 4, we 
solve for ^ as 

= ^^^x + - J (^^x + ^j,?/) dt, 
= x<i/^ + y^y + f{s). 

Here we used + = 0. At t = 0, ^(1, s) = log(l + s) + /(s) and 



vl>,(l,s) 



9^(1, s) 



+ fis) = 0. 



5s s + 1 

This implies that /(s) = — log(s + 1) + constant. We will later show that 
must vanish at x = 1, y = (i.e., t = 0, s = 0) so that this constant = 0. 



Therefore, we have (|3.30|) . 

We plot the rays using (|3.33|) and (|3.34|) in Fig. 4 and Fig. 5. In Fig. 
4 we plot the rays for s > and t in the range < t < tmax = l/"^- 
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Att — we have x — 1 and at t — tmax{s) the rays reach x — 0. We 
note that develops a singularity at t = tmax{s). We continue the rays 
for t > tjnax in Fig. 5. There is a cusp when x — 0, which happens at 
t — tc{s) = 1 + 1/s > tmax{s). Att — tc, X + y — and thus the locus of the 
cusps is on the line y = —x. By comparing Fig. 2 and Fig. 4 we see that the 
main difference in the rays occurs in their behavior near y = 0, and indeed 
the asymptotics will be much different for y ^ and p t 1/ than what we 
found for w and p < 1 in section 4. Near y = the ray expression 
breaks down and we will analyze these cases in subsections 5.2 - 5.5. 



2.0 



1.5 



1.0 



0.5 








Fig. 4 A sketch of the rays from x — 1. 
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-1.0 



-0.5 



0.5 



1.0 



Fig. 5 A sketch of the extended rays from x = 1. 



We next calculate IC{x,y) in (|5.4|) . In the same manner as in section 4 
we obtain the following 'transport' equation: 



{ 



- a] {x + y) + 1 



e-"^ + {^^xx + l) e*^ + [l%y + l) e*^ 



-2 + a{x + y)})C= [{x + 2/)e-*- - xe^^] IC^ - ye^yJCy. (5.23) 



The right hand side is xK-x + y^y = and thus, from (|5.8|) , (|5.9|) , (|5.11|) , 
(|5ia , and d^ISl) - (l430l) . we can write (ISlISl) as 



i-^ + l(vl/, + + a{x + y)il- e-*^) + ^e-*^ 



(5.24) 



Integrating (|5.24|) using (|3.31|) for we obtain 



]C = ICois)\j\-'/'exp i ^^i±ly ]\s-st + ll^/^exp 



ast 



2(s + 1) ^ 2 



(5.25) 
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where /Co(s) is a function of s. Using (|3.33|) and (|3.34[) , we obtain j as 

s{s -st + 1) 



{2s{s + 2) - [s{s + 2)t^ -2t + s{s + 2) - l]e~*} . 



(5.26) 



Therefore, 



/Co(s)(s+l)3/2(s-st + l)e 



t/2 



v/s(l - st)(se* - St + l)y2s(s + 2) - + 2)t^ -2t + s{s + 2) - l]e^ 



X exp 



2{s + l) 

We determine /CqI-s) by matching V{x, y) to the ray solution in section 
4. Setting p = 1 — a6 in (|3.2[) and expanding the result for 5 — > yields 



(5.27) 



— a5 , 

p=i \dp 



p=i 



p=i 



p=l 



- log(s + 1) - a6{(j)sSp + 4>ttp + 



(5.28) 



We can easily check that when p = 1, (p^ = (f)y = "^y, and (|3.5|) and (|3.6|) 
reduce to (|3.33|) and (|3.34|) respectively. After some calculation we find that 
when p = I, 



^sSp + (t>ttp + 



,Sp + (j)ttp + Xp(f)^ + x(j)^p + yp(j)y + y(f)yp + 1 



2{s + l) 



- st + s + l. 



(5.29) 



Thus (|5.28|) becomes, as 5 ^ 0, 

~ + y^!y - log(s + 1) + a(5 
and expanding (|3.7[) as 5 — gives 



St 



2(s + 1) 

g3/2 (g _ + 1)3/2 

s + iy Vse* - St + 1 ■ 
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s + l) 



(5.30) 



(5.31) 



Therefore, from (|3.30|l , (|5.27|) , (|5.30|l , and (|5.31|) , the matching suggests that 



as 



3/2 



and thus we obtain (|3.35)) . 



The ray solution V{x, y) ~ 6IC{x, y) exp[\l/(x, y)/S] is valid in the inte- 
rior of the state space and also near x = 1. However, the expansion breaks 
down near the boundaries x = and y = 0, and near the comer points 
(x, y) = (1, 0) and (0, 0). We analyze these regions separately in the follow- 
ing subsections. 

Near x = 1 we can simplify the ray solution to 

Say 

7r(fc,r) = V{x,y) ~ —^exp[m{x - l)log(?/ + 1) -a(?/ + 1)], (5.32) 

so that the exponent is 0(1) for k = m — 0(1) (i.e., l — x = 0(m"^) = 0{6)). 
The total probability mass in this region is given by 

M+ = Prob[Ni = m-0{l), N2 = 0{m)] 

poo o^, 7, 

~a/ J2(y + l)--^e-'^^y'''^dy = e-\ 

In subsection 5.2 we shall consider the scale < x < 1 and r = 0(1), and 
obtain the mass in this range as 

M- = Prob [0<Ni<m, iVa = 0] ~ 1 - e'^, 

which complements Ai+. Thus for m ^ oo and in the heavy traffic limit, 
there will either be no secondary spaces occupied, or the number of occu- 
pied primary spaces will be close to the maximum m and then there will 
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be a large number (0(m)) of secondary spaces occupied. From (|5.32|) we 
can also infer the conditional limit laws 



Prob[Ni = m - n\N2 = my] ~ y{y + 1) " \ y > 0, 

P h\Ar \Ar 1 y(l/+l)'""^e-"^ 
i^roo [J\2 = ^7i?/|A/i = m — ~ , 

^71 



with 



C 



-du. 



'0 

This may be written in terms of incomplete Gamma functions as 



r(l — n,a) — r{—n, a) 



5.1 Boundary layer near x = 

We consider the scale k = 0(1) and y > 0, which corresponds to x = 0{5) 
and r = 0(5^^) = 0{m). From (|3.31|) and (|3.33|) we see that has a 
logarithmic singularity as x 0. We set 



n{k,r) = ^exp 



(5.33) 



and rewrite (|2.3|) as 



(2-a5)[y+(A;+l)(5]e*(°'^)/^5,(y; S) = S{l-aS)[y+{k+me''^'''^^'Sk-i{y; S) 
+ ik + l)e*(°'^)/^5fc+i(i/; S) + {y + 5)e^^''y^'^"Sk{y + 5; 5), 

A; > 0, y > 0. (5.34) 
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The boundary condition at = can be written as 



(5.35) 

We expand Skiy; 6) in the form 



S,{y;S) = S,{y) + 6S;^\y) + 0{6'). 



and let 5 ^ in (|5.34|) and (|5.35|) to obtain the following equation for the 
leading term Sk{y) 

[2 - e*-(0'^)] ySkiy) = {k + l)Sk+i{y), k > 0, y > 0. (5.36) 

The general solution to (|5.36|) is 



and thus 



(5 



^(0,1/) 



(5.37) 



In the same way as in section 4 we determine the constant u and the func- 
tion So{y) by matching (|5.37[) asymptotically to the ray solution V{x, y) ~ 
5IC{x, ?/)e*'^^'^)/^ in an intermediate limit (x ^ but k = x/5 ^ oo). Using 
Stirling's formula for k\ we can write the expansion in (|5.37|) for /c ^ oo in 
terms of x, as 



7r(A;, r) ~ ' exp 



TCX 



i 1^(0, + X + xlog (^^) + xlog [2 - 6*"^°'^)] } 

(5.38) 
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We expand the ray solution V{x, as x — > 0, inverting the transfor- 
mation from (x, y) to ray coordinates along x = 0. When x = we have 
t = tmaxis) = l/s. Then from (|3.34[) we find that s and y are related by 



y = I (s + 1). We thus define s = s{y) by 



s{y) = 2 



(5.39) 



Note that this is the same as (|3.13|) with p = 1. Thus, a ray that starts from 
(x, = (1, s) hits the y— axis at the point (0, y). Using (|3.30|) and (|3.32[) 



and evaluating these expressions at t = tmax and s = s we obtain explicit 
expressions for \E'(0, y) and ^y(0, y): 



vl/(0,y) 



s + 1 



-tr 



log(s + l), 

1 

s 



which can be rewritten as (|3.38|) and (|3.39|) respectively. Also, from (|3.31|) 
and (|3.33|) we obtain 

;i + s-st)(l + 2s-st-se-* 



log X + log 



s + 1 



(5.40) 



which is exact for all (x, y) and indicates the logarithmic singularity in 
as X — > 0. By expanding (|5.40|) for s ^ s and t t^ax we obtain 



r fi2 



which integrates to 



log X + log 



+ o{l] 



s + 1 

logx + logy + log [2-e*^(°'^)] +o(l), 



(5.41) 



^(x,?/) = ^(0,2/) - xlogx + x + xlogy + xlog [2 - e'^^^^'^)] +o(x). (5.42) 
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The above agrees precisely with the exponential terms in (|5.38[) and thus 
the expansions will asymptotically match provided that 0=1/2 and iSo (?/) = 
lim^j^o [VStts /C(s, y)\ . By expanding j = XtUs — XsVt and then /C in p.35[) 
for s ^ s and t — > tmax, we obtain 



as 



(s + l)Vl+2 



exp 



a s + 



2(s+ !);_ 



(5.43) 



which can also be written as (|3.40[) in terms of y. We have thus determined 
V and iSo(|/) in (|5.37|) , and verified the matching between the {x,y) and 
(A;, y) scales. 



5.2 Boundary layer near y = 

We consider small values of y and use the original discrete variable r 
y/5. We set, for < x < 1, 



(5.44) 



and 



7r(A;,0) = (l-p)p'= + m'^iv^e-^^(^)Qo(a;;5), r = 0, 



(5.45) 



where 7i(s) is a positive function of x. The forms of the expansions in 
(|5.44|) and (|5.45|) are indicated by the behavior of the ray expansion asy ^ 
0, for < X < 1, which we shall discuss shortly. We also note that in the 
present heavy traffic limit 



(5.46) 
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With BA^ and (|5^ , from dZS]), we obtain 



(2-a5)[x+5(r+l)]e-v^^(^')Q,(x; 5) = {l-a6)[x+6{r+l)]e~^'^^''~^^ Qr{x- 



(5.47) 

and 

(2 - a5){x + 5)e-^^(")Qo(a;; 5) = (1 - a5)(a; + 6)e-^'^^''-^^ Qo{x - 6; 6) 
+ {x + 5)e-^'^^''^^^Qo{x + 5;5) + 5e-^'^^^^Qi{x;5), r = 0. (5.48) 

Expanding Qr{x; 5) and Qo{x; 5) as 

Qr{x;S) = Qr{x) + SQ^\x) + OiS^), r > 0, 



we obtain from (|5.47|) at the first two orders {0{5) and 0(5^/^)) the follow- 
ing equations 

[n\x)r 

2n'{x)Q'^{x) 



1 

X ' 



-n'(x) -aH'ix) -n"(x) 



L X 



Qrix) 



(5.49) 
(5.50) 



Taking 7-^'(x) < 0, from (|5.49|) we obtain 7Y(x) = —2y/x + constant. By 
matching to the ray solution at x = 1 (corresponding to t = 0), where 
\E' = 0, we must have H{1) = so that 



n{x) = 2 (1 - v^) . 



(5.51) 



Using (|5.51|) we solve (|5.50|) to obtain 



(5.52) 
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Here Q*(r) is a function of r. Thus, for r > 0, the boundary solution near 

?/ = is 

We determine (r) and vi by matching this boundary layer solution 
to the ray solution V{x, y). We expand the ray solution V{x, y) as y 0. 
We let s ^ 0, t ^ oo with st fixed and < st < 1. Then (13301) . (13331) , 
(I334b , and (l335l) give 

(1 - 2s(l - St) 



X ~ 

2/ ~ 

/C ~ 

V]/ ^ 



s + 1 

^2 



+ 



S + 1 



S + 1 



- Vl — st exp 



—a 1 — st + 



;i-st)^iog 1 + 



1 - st 



s log 



1-st 



se' 



log(l + s) 



s2[log(l-st)-logs-2t]. 



(5.54) 
(5.55) 
(5.56) 

(5.57) 



We can invert (|5.54[) and (|5.55[) and write s and t in terms of x and for 
?/ ^ 0, as 



1-st 



X 



1 — ^/x 1 



(5.58) 
(5.59) 

(5.60) 
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Using (15.581) - (I5.6UI) in (15.561) and (15.571) we then obtain 

51C ~ ^a;i/4e-(-+i)/2, (5.6I) 

vi> ^ _2^(i_y^) + |iog(^^^, (5.62) 

and thus near y = the ray expansion becomes 



Therefore, comparing (|5.53|) to (|5.63|) we conclude that Q*(r) = |e "/^r ''/^ 



and ui = —1. We have thus obtained (|3.41[) . 
For r = 0, we rewrite (|5.45|) as 



-2y?S(l-v/5) 

7r(fc, 0) ~ (1 - p)p^ + ^ Qo(x). (5.64) 



Then, from (|5.48|) , we obtain to the leading order 



Q,(x) = = -Q.(l)a;3/^e-'^^/2 = -^x'/'e-^^^^'^/\ (5.65) 



and thus obtain (|3.42|) . 



5.3 Corner layer near (x, y) = (0, 0) 

We consider {x, y) near (0, 0) and go back to the original discrete variables 
(/c, r), with k,r = 0(1). We also set 

7r(A;,r) = m-^/^e-^^ r(A;, r; m), r > 0, (5.66) 
vr(A;,0) = (1 - p)p'^ + m'^/V^^ T(fc, 0; m), r = 0, (5.67) 
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and note that T(0, 0; m) = 0. The scaling in (|5.66[) and (|5.67|) can be in- 
ferred by expanding (|3.41|) and (|3.42[) as x = 6k = k/m ^ 0. Then the 



balance equation (|2.3[) becomes 

(2-a5)(A;+r+l)e"^^r(fc,r; m) = (l-a(5)(A;+r+l)e-2^ T(A;-1, r; m) 

+ {k + l)e-^^ r{k + 1, r; m) + (r + 1)6" V"^('-+i) r(A;, r + 1; m), 

(5.68) 

for A; > 0, r > 0, and 

(2 - a5)(fc + 1) T{k, 0; m) = (1 - a(5)(fc + 1) T{k - 1, 0; m) 

+ (A; + 1) T(A; + 1,0; m) +T(A;,1; m), A; > 0, r = 0. (5.69) 

The boundary condition at A; = is 

(2-a6) e-2v^T(0,r; m) = e-^"^ T{l,r; m)+e- V™('^+i) T(0,r+1; 

r + 1 

(5.70) 

for r > 0, and the comer condition at (k, r) = (0, 0) is 

= T(l, 0; m) + T(0, 1; m). (5.71) 

Expanding T(k,r; m) as T{k,r; m) = T{k,r) + o(l) we obtain from 
(|5.68|) - (|5.70|) the following equations for the leading term 

2(A; + r + 1) T{k, r) = {k + r + 1) T{k -l,r) + {k + 1) T{k + 1, r), 

A; > 0, r > 0, (5.72) 



2(A; + 1) T{k, 0) = {k + 1) r(A; - 1, 0) + (A; + 1) r(A; + 1, 0) + T{k, 1), 

A; > 0, r = 0, (5.73) 
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2(r + 1) r(0,r) = r(l,r), r > 0. (5.74) 
We shall first analyze (|5.72[) and (|5.74|) for r > 0, and then solve (|5.73|) for 

r(fc,o). 

Forr > we use the generating function ^ (z) = Q{z]r) = XlfcLo ■^(^' ^)'^'^- 
Multiplying (|5.72|) by 2;'^ and summing over k > gives 



^ - lyg'iz) = [-{r + 2)z + 2(r + 



(5.75) 



The general solution to (|5.75|) is 



where ^*(r) is a function of r. Inverting the generating function, we can 
write 



r{k,r)=g,{r)- 



1 /• z-^-^ 



f/^^-^'^dz, r>0 



(5.76) 



27ri / (1 - zy+^ 

where the integral is a complex contour integral along a small loop around 
z = 0. We obtain ^=„(r) by asymptotic matching. Letting x ^ {x = k/m) 
in (|3.4H) , and invoking the asymptotic matching condition shows that 



We now expand (|5.76|) as A; ^ 00. This will verify that (|5.77[) is satisfied and 
also determine the function ^*(r). As A; 00 the behavior of the integral 
in (|5.76|) will be determined by a saddle point near z = 1. We write the 
integrand in (|5.76|) as 

exp[H{z)] _ 1 



z(l - zY+^ 



'l-zY+^ 



exp 



—klogz + 



1-z 



(5.78) 
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The saddle point(s) satisfy H'{z) = so that 



k 



(5.79) 



z (l-z)^' 

(|5.79|) defines the saddle i as a function oik/r, and it is given explicitly by 



r ^Jr{Ak + r) [r 

'^2k — 2k — '-y^'^ 



oo. 



(5.80) 



Evaluating the integral along the steepest descent path through the saddle 
point we obtain, as A; — > oo, 

r(k, r) ~ -J—r--/^-y^k'/^+'/''e'/^e^^g,{r). (5.81) 



Therefore, BIT} and BM\ suggest that ^*(r) = ay^ ^3/4g-(a+r)/2^ 
thus, for r > 1, we obtain (|3.44|) . 



Next we show that (|3.44|) , when expanded for r ^ oo asymptotically 
matches to the expansion in the boundary solution near x = 0, as y ^ 0. 



Letting y^Oin (l338l) - (l3^ we obtain 



+ 0(1). 



Thus as (?/ = r/m) the asymptotic behavior of (|5.37[) is 



k\ 



(5.82) 



This is an approximation to 7r{k,r) that applies in the matching region 
where A; = 0(1), r — > oo, r/m — > 0. 
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Next we expand the corner approximation (|3.44[) as r ^ oo. Expanding 
the integrand in (|3.44|) around z = by setting z = u/r and evaluating the 
resulting integral, using 



u 



-k-12u 



shows that (|3.44|) , as r — cxo with k = 0(1), agrees with (|5.82|) . 

We verify that the comer approximation (|3.44|) matches to the ray so- 
lution 5}C{x, y) exp[m\E'(x, y)], in an intermediate limit where k,r ^ oo but 
X = k/m, y = r/m ^ 0. We expand (|3.44|) for k and r simultaneously 
large, writing the integrand as 

exp[H{z)] _ 1 



z{l-zy z{l-z)' 



exp 



—k log z — r log(l — z) 



and again using the saddle point method. The saddle point(s) now satisfy 
H'(z) = so that 



k 
r 



1 + 



1 



(5.83) 



l-z \ 1-z^ 

(|5.83|) defines the saddle 5 as a function oik/r, and it is given explicitly by 



z = l- 



k + r 



We note that z as k/r and z 1 as k/r oo. Evaluating 
the integral along the steepest descent path through the saddle point we 
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obtain 



7,{k,r) ~ « ^-5/4^1/4g-(a+r)/2g-2v^g./(l-.-)^^(L^ 



- ^-5/4^1/4g-(a+r)/2g-2v^g^r(fe+r) 

-fe-V2 /^^^\r-/2+l/4 



This approximation applies in the matching region where k,r ^ oo, but 

k/m,r /m —>■ 0. 

In the ray solution we let x ^ and y ^ 0, which corresponds to s ^ 
and t — > oo, with x/y fixed (corresponding to 1 — st ~ s). We find that in 
this limit 



1 y/x + y 



1 - St ~ V^ + y - Vy, s y/y, tr^ — —— + 1, 

and thus, from (l330ll - (l332l) and (|335)) . obtain 





~ log ( 




~ log ( 




~ X loe 



y/xTy 



y/x + y-y/y 

VxTy\ 1 _^ yx + y _ ^ 



Vy J Vy Vy 

VxTy \ r-, — ■ — r y , ^ ( VxTy 
r Q ,-a/2 + 1/ 

Therefore, as (x, y) (0, 0), in terms of a; = k/m and y = r/m, 



V{x,y) ~ ^ ^-lg-(a+r)/2g-2V^gy;(fcTO 
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We can easily check that (|5.85[) agrees precisely with (|5.84[) . 

It remains to compute T{k, 0), by solving (|5.73|) and (|5.71|| . We use the 
generating function 

oo 
k=0 

multiply (|5.73)) by z'', and sum over A; > 0, to obtain 

{z - + 2{z - l)go{z) + g,{l){l - 2)-3ei/(^-^) = 0. (5.86) 



The solution to (|5.86|) is 

,1/(1-^) 



ze 



z-l) 



-{a+l)/2 



z 1 - z)-=e 



3^1/(1-^) 



and thus. 



mo) = — (1^ z-'-'go{z)dz 



— ai/vre 



-(a+l)/2_ 



(5.87) 



e^/^'^-'Uz. (5.88) 



2m J (1-2)3 

To verify the matching between (|5.67|) and (|3.42|) we evaluate the inte- 
gral in (|5.88[) as A; ^ oo. We rewrite the integrand in (|5.88[) as 



exp[Ho{z)] 



1 



exp 



—klogz + 



1-z 



(1-Z)3 (1-^)3 

and use the saddle point method. The saddle point(s) zq satisfy Hq{z) = 
so that k = z/{l — zY and hence 

. ^ V^k + T 1^1 
= ' ^^Tk^'-7=k 
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Evaluating the integral along the steepest descent path through the saddle 
point and using the result in (|5.67|) we obtain 



TT 



{k, 0) ~ (1 - p)p'= - -e-"/2m-5/4^3/4g-2v/Hg2v^_ 

2 



(5.89) 



As X ^ {x = k/m) the asymptotic behavior of (|3.42|) agrees with (|5.89|) . 
This completes the analysis of the corner range where {x, y) ^ (0, 0). 

5.4 Analysis near the comer (x, y) = (1, 0) : 

We examine the problem for x = 1 — 0(^/6) and y = 0(^/6), which corre- 
sponds to k = m — 0{^/m) and r = 0{^/m). We let 



(5.90) 



and set 

7r(fc,r) = 5). (5.91) 

The scaling in (|5.91|) can be inferred by expanding (|5.32[) as x = 1 — 1 
and y = R\f5 — > 0. Then, for < ,^ < 1 and i? > 0, (|2.3|) becomes 



(2 - a(5)r](^, R; 5) = - V6, R; 6) + {1 - a6)n{^ + V6, R; 5) 
{RV6 + 6)n{^,R + V6; 6) _ rV6 - V6, R; 6) 
l + {R-0^ + 6 l + {R-0VS + 6 



. (5.92) 



The boundary condition (|2.7[) becomes 

1 + 6 



1 + RV5 + 5 



n{-VS,R; S) = (1-6) n{0,R-VS; S), R>0. (5.93) 
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Expanding R; 5) as 

n{C,R; s) = n{^,R) + o{Vs), 

from (|5.92|) and (|5.93|) we obtain to the leading order 



+ % + + 1]^) = ; ^, R>0 (5.94) 
i? + fig - = ; ^ = 0, i? > 0. (5.95) 

We must thus solve a parabolic PDE in the quarter plane, subject to 
an oblique derivative boundary condition along ^ = 0. While we were 
not able to solve this problem exactly, we shall use asymptotic matching 
to infer various properties of R) as ^ and/ or R become(s) large, or if 
R ^ 0. We shall also obtain an integral equation for the boundary val- 
ues f2(0, R), and show that the problem can be reduced to the heat equa- 
tion with a moving boundary. We comment that the PDE in (|5.94|) is not 



separable due to the term Ril^. When considering the analogous infinite 
server model (see [8J) Knessl obtained, in a certain heavy traffic limit, a 
problem very similar to (|5.94|) and (|5.95|| . However, there the term RQ^ 



was replaced by ^fi^, so that the PDE was separable. Then despite the 
oblique derivative boundary condition Knessl could solve this problem 
explicitly, in terms of contour integrals of parabolic cylinder functions. For 
the infinite server model the corresponding probability 7r(k, r) on the scale 
(|5.90|) was 0(m^^) and most of the probability mass accumulated on this 



scale. In the present processor sharing model we have, in view of (|5.91|) , 

68 



7i{k,r) = 0(m~^/^) for k = m — 0{^/m) and r = 0{^/m). Thus the total 
mass is roughly 0(m^^/^) x 0{^/m) x 0{^/m) = 0(m^^/^). Indeed we al- 
ready showed that for the PS model with m — > oo and 1— p = 0(m^^), most 
mass occurs either along r = 0, k = 0{m) or r = 0{m), k = m — 0(1). The 
PS model thus, on the R) scale, leads to a more difficult mathematical 
problem, but one whose analysis is perhaps less critical, due to the small 
probability mass. 



We first examine (|5.94|) and (|5.95|) for ^ and/ or _R oo. This will estab- 



lish the asymptotic matching between the R) scale and the ray expan- 
sion on the (x, y) scale. We introduce the small parameter e > 0, setting 

e e 

and 

n = £''*e-^(«''^*)/^'i3(r, R*) [1 + 0{e^)\ (5.96) 
From (|5.94|| and (|5.95|) at the first two orders (0(6^^) and 0(1)) we obtain 



AI, + R*{A^*+Ar*) = Q, (5.97) 

2^5* iSg. + R*{B^* + i3jj.) + (^5*g. + l)i3 = (5.98) 

and when ^ = (ESSj) leads to 

- An, + i?* = 0, (5.99) 

i?^. - Br* = 0. (5.100) 

We are thus using a ray expansion on the PDE (|5.94|) . We again use the 
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method of characteristics, writing (|5.97|) as 
where 

J^* = Al. + R* (^e + ^fi* ) • (5-101) 
The characteristic equations are 



duj dA(^* 
dR* dT* 



duj dAfi* 



2^5* - R\ (5.102) 
-R\ (5.103) 



doj 

dA^* _ dr 



dcu d^* 



2Al, - R* (^e + Ar*) = -Aj, , (5.104) 
0, (5.105) 

-Ai:,-AR*. (5.106) 



duj dR* 

Letting R* = u at u = (using rays start from {R*, ^*) = {u, 0) at = 0) 
from (l5J02)l - (|5l06)l and (15:99)) we obtain 



R* = eR = ue-'^, (5.107) 
f = £^ = u [2uj + e-'^ - 1] , (5.108) 
A = -u^uj. (5.109) 

Here we also used ^(0, 0) = to determine an integration constant. We 
rewrite (|5.98|) as 

^ = (^ee + 1)^ (5.110) 
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The general solution to (|5.100|l and (|5.110|) is 



13 = Bo- 



ue 



(5.111) 



V2a7 + T 
where Bq is a constant. 

To determine Bq we examine the asymptotic behavior of the solution 
(|5.96|) as ^ ^ and R ^ oo, which corresponds to m — > cx) and a; — >^ 0. 
From (I5l07b - (I5J091) and (|5lTT1) we obtain 



Thus, from (ISlgTl) and (15:961) 



7r(fc,r)~53/2^^*+%i?e-^. 



(5.112) 



Letting |/ — »• (|/ = v^i?) and x = 1 — in (|5.32|) gives 

7r(A;,r) ~53/2^i?e-"e-«^, 
and matching thus forces z/* = — 1 and Bq = ae~°. Therefore, 



(5.113) 



= exp (- 



-2 2 \ 
-€ U to] 



Dp 2a; 

ae " , = exp i—R cue 



(5.114) 
(5.115) 



Here we used u = sRe"^. This gives the behavior of i?) as i? oo 
withO < ^/R < oo. 

We can rewrite (|5.115|) in terms of the Lambert 14^-function. From (|5.107|) 
and (I5J081) we obtain ^/R = {2uj + - 1) = e'^(2c<; - 1) + 1, so that 



|-l = 2v/iL-l)e--V2 



(5.116) 
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Setting W = uj-\l2 (l5J16ll becomes 



We 



w 



so that u = 1/2 + W, where W is the Lambert function. Thus (|5.115[) 
becomes 



4R^2{W + 1)W- 



exp 



21 W + 1/2 



(5.117) 



We next show that (|5.117|) or (|5.115|) also follows by expanding the ray 
expansion on the {x,y) scale, as {x,y) (0,0) along lines where y/x is 
fixed. We expand the ray solution V{x,y) as s ^ with t fixed. Letting 
s ^ with t fixed in (|3.33|) and (|3.34[) gives 



x[s 



t) = 1 + s (1 - 2t - e~*) + 0( 



y{s,t) = se~* + 0(s^ 



and thus ye* ~ s and 

1 — X — y 1 — X 



s ~ 



2t - 1 ' y ' " ~ \ 2 
From SSJTE^ . (|330b - (l332b , and SSlTh we obtain 



(5.118) 



X 



s + sUt 



-s + 



~ > --s]il-x-y) + s^[ix + y)t-l] 



(1 - X - yf 
(2t-l)2 



t + (l-x-y)t 



1 — X — y 



(5.119) 
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and 



ae ^ {1 — X — y)" 



(5.120) 



Since {1 — x)/y = ^/R, in view of (|5.116|) and (|5.118|) we can identify t, as 
t ~ u; in this limit. Therefore, letting x = 1 — and y = \f5R in (|5.119|) 
and (|5.120|) we find that V{x, y) agrees with (|5.115|) . 

We next examine the asymptotic behavior of (|5.115|) as i? ^ and — *• 
oo (corresponding to m ^ and uj ^ oo with uuj fixed). From (|5.107|) and 
(|5.108|) we obtain 



u 



2a;' 



e 1 /^^ 1 

, uj ~ log — . 

2ujR' ^\R 



Then (|5.114|) becomes 



R 



logR\ ^''^exp 



41og(e/i?)J 



(5.121) 



If we fix R and let ^ ^ oo (corresponding to u,uj oo), similarly we 
obtain 



■5/2 



exp 



(5.122) 



4iog(e/i?). 

A more uniform result can be obtained by replacing log^ in (|5.122|) by 
log(^/i?), which will contain also (|5.121|) as a special case. 

We next examine (|5.94|) and (|5.95|) as ^ i 0, R ^ oo with ^R fixed. We 



shall thus obtain some of the higher order terms in (|5.113|) . We set 

R* 



R 
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and expand R) as 



(C, R*) + en'-^^ ic, R*) + e^n^^^ (C, R*) + 0{e') 



Then, from (|5.94[) we obtain for the first three orders {0{e 0(1), 0(£:^)) 



)(0) 



^^o^o H~ R ^ 



(0) 

(1) 



(5.123) 
(5.124) 
(5.125) 



and from (|5.95|l we obtain the following boundary conditions at = 

0, 



nfj + R*Q^^^ 



(5.126) 
(5.127) 
(5.128) 



Solving (|5.123|) - (|5.128|) recursively we obtain after some calculation 













2 








"i?*(^^)4 






8 


6 



1 



(5.129) 
(5.130) 

. (5.131) 



2R* (i?*)2 (i?*)3_ 
We note that (l5J23ll and (l5l26l) imply that = e-«*^7o(^*), and we 



ultimately determine /o from the solvability condition for (|5.125|) . Then to 
obtain Q'^^^ we must also consider the equation for etc. Therefore, as 
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n and i? ^ oo with = 0(1), 



' i?3 i?2 2i? 6 8 



0{R-^)]. (5.132) 



Next we analyze (|5.94|) and (|5.95|) using a Laplace transform. We set 
n{^, R) = V{^, R)/R. Then we rewrite dSlM)) and (15:951) as 



+ i?(Pg + Vr) = ; ^ > 0, 



(5.133) 
(5.134) 



We assume that ^ as i? ^ and use a double Laplace transform, with 

POO 

U*{R;a)= V{i,R)e-''^ di, (5.135) 
Jo 



and 



roo 

U{a,(3) = / U*{R,a)e-^^ dR, 
Jo 

POO POO 

= V{^,R)e-"^e-'^^ d^ dR. 

Jo Jo 



(5.136) 



Taking the Laplace transform of (|5.133|) over ^ gives 



a^U*{R; a) + aRU*{R] a) + RU*j^{R] a) = aV{0, R) + ©^(O, R) + RV{0, R). 

(5.137) 



Taking the Laplace transform of (|5.137|) over R we obtain 

a'U - (a + f3)U, = e"/^^ |aP(0, R) + r\-^ i^^)] } 



dR. 
(5.138) 
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Here we also used (|5.134[) to eliminate V^{0, R). Integrating by parts in the 
right-hand side and dividing both sides of (|5.138|) hy —{a + (3) we obtain 



U 



R{a + (3) 



l)V{0,R)dR, (5.139) 



a + P y JQ 

and we can rewrite (|5.139|) as 

POO 

L h Jo 

We assume that the real part of is negative, i.e. 3?(q:^) < 0, integrate 
(|5.140|) over {—a, (3), and obtain 



^-a-/3] V{0, R)e-^^{a + (3y^ dR. (5.140) 



U{a,f3) = 



1 r f a + l3 

a + (3 7_Q 7o V " + ^ 



— -a-z] ©(0, R)e-^^ dR dz 
R 



1 



1 + 9 



--{a + p){l + e) 
X 



P(0, ;^)e-«e-("+^)^^ dx dd. 

(5.141) 

Here we changed variables with z = {a + 13)9 + (3 and R = x- Inverting 
(|5.141|) over (3 gives 



U*{R;a) 



-1 



1 + ^^ 



X 



1 aO \^f^ R 



i + e 



1 + ^ 







_ _ _ p 0, ^ - 0, ^ 



X X/ X 



de 
de 



dx, 

(5.142) 



where we set = — 1 and 9 = R/x- 
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Inverting (|5.142|) over a we obtain 



l^i^.R) = f {7^ f exp[a(x-i? + + «'logX-«'logi?] 

Jr 27r2 J Br 

X [(1 - ax)V{0, x) - X'D^iO, x)] da} ^ 

ix~R + 0'' 



41og(x/i?) 



1 ^ x-R + ^ \ P(o,x) 
X 2\og{x/R)J xv^MxAR) xV^og{x/R) 



dx 



ix-R + O' p(o,x) 

20Fy« 4[\og{x/RW/' """^ 



1 d [ l?(0,x) 

'^Vt^ Jr dx [ X 



41og(x/i?) 
(x-^ + Oni dx 



(5.143) 



41og(x/i?) 



(5.144) 



Integrating the second integral by parts, we obtain an alternate expression 
for V{^, R) as 



'D{^,R) 



'^V^.Jr 



ix-R + O' 



'D{0,x) 



A[log{x/RW/' 2[log(x/i?)]3/2j X' 

ix-R + 0' 

41og(x/i?) 



X exp 



dx. (5.145) 



Next we derive an integral equation for V{0, R) by letting ,^ ^ in 
(15.1431) . We can rewrite (I5.143D as (13.471) . If ^ = 0, the second integral 
has a non-integrable singularity at the boundary x = R> due to the factor 
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[log(x/-R)] If we set x = r the second integral becomes 



exp 



2^3/2 

1 f°°V{0,Re' 



Av 



dv 



4 A io 



w 



3/2 



exp 



R^je^^' - l + i/Rf 



dw (5.146) 



where f = wS, . Then letting ^ | in (|5.146|) we obtain 



40F 



V{0,R) 
20F 



«-V^e-^ci« = (5.147) 



Therefore, from (13^ and (I5l47|) we obtain (13^ . Since ^(0, i?) = /?fi(0, R) 
we thus also have the following integral equation for 0(0, R) 

ix - R? 



Rn{o,R) = -= 



v^Jr ^J\og{x/R) 



exp 



x-R 



41og(x/i?)J 



r](o,x)-f^x(o,x) 



dX- (5.148) 



_21og(x/i?) 

We can also analyze the integral equation (|5.148|) asymptotically. Con- 
sider the limit R oo. Then assuming that f2(0, R) has mild (e.g. alge- 
braic) growth, the function exp{ — (x — -R)^/[4 log(x/-R)]}fi(0, x) becomes 
sharply concentrated at x = ^/ which is the lower limit on the integral in 
(I5ll8ll . We thus set x = ^ + V/R in (I5ll8ll and note that 



ix-Rf V' 



+ 0{R~ 



(5.149) 



log(x/i?) 2i?2 i2i?4 

Then if we evaluate (|5.148|) by an implicit form of Watson's lemma, to lead- 
ing order the right-hand side becomes i?fi(0, R), which is the same as the 
left-hand side, but this argument does not determine fi(0, R). However, by 
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considering higher order terms in the expansion of the integral, we find af- 
ter some calculation that the right-hand side of (|5.148|) has the form, with 

n{o,R) = no{R), 

+ [SR^n'oiR) - SGR^n'oiR) + 88RnQ{R) - lOOQoiR)] + 0{R-'^no). 

(5.150) 

2n" 



We note that if VIq has algebraic behavior as i? ^ oo then Hq, RVIq, R^Q, 



0; 



all have the same order of magnitude. The expansion in (|5.150|) is in pow- 
ers of R^^, but the coefficient of i?^fio turned out to be zero. Setting 
(|5.150|) equal to RiloiR) we see that asymptotically the integral equation 
BTim can be approximated by 2R^n'^{R) - 4:Rn'^{R) + 4fio(^) = 0. This 
is an ODE of Cauchy-Euler type that admits the solutions fio(-R) = R and 
^o{R) = R^. But by asymptotic matching to the ray expansion (x = 1, 
y > 0) we know that ^oiR) ~ ae~'^R, which precludes the second solution. 



Then from (I5l50l) we can conclude that if n^iR) = ae~''[R+cR-^ + 0{R-^)] 



as R oo, then c = 1. This result agrees precisely with (|5.132|) with = 0. 



The third (— 6-R ^) term in the expansion of flo{R) can also be obtained by 
explicitly evaluating the 0{R^'^ilo) terms in (|5.150|) . 



Now consider the limit ^,R oo with ^/R fixed, and equation (|3.47|) 
We evaluate the integral by an implicit form of the Laplace method. Scal- 
ing X = Rv the integrand will be maximal where 



d f [R(v-l)+^f 
dv \ \ogv 
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0, 



or 



-I- 



2v log V = 

Letting v = -\/ee" this equation becomes equivalent to (|5.116|) , with u = 
W{{^/R — l)/{2y/e)) being the Lambert-iy function. Then if we evaluate 
(|3.47|) by the Laplace method we obtain, after multiplying by R, 



27r 



RQiC,R) ~ exp [-R^w.vl] J — 



X 



[R{v, - 1) + ^]Rno{Rv,] 



d 
dv 



(5.151) 



where f* is the location of the maximum, w^: = logf* and = H"{v^) 
where H{y) = [R(v — 1)+^]'^/ (4 log v) . We see that the first term in the right- 
hand side of (|5.151|) dominates the second, and after some calculation we 



find that * = R^{1 + 2w^)/{2w^) so that (l5l5T1) becomes 



n{^,R) ~ QoiRv. 



exp [—R^w^vl) . 



(5.152) 



~ ae~"-Rv^:, so that (|5.152|) agrees precisely with 



Now by matching ^7o {Rv* 
(|5.115|l , after we identify w^, c<j, f ^, ^ e'^. 

Thus (|5.145|) or (|3.47|) yields quite a bit of information as i? ^ oo, since 
we can then localize the integral operator. However, the opposite limit 
i? ^ 0+ seems much more difficult as then the global nature of the op- 
erator persists, and the asymptotic evaluation of the integral seems to re- 
quire considering separately the contributions from different ranges (such 
as X ~ -R/ X = 0(1) and x = 0{R)). If x = 0(1) the contribution will 
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of course involve 1^(0, x) (or ^o(x)) for x = 0(1), which we do not have 
explicitly. 

We can also derive a heat equation from (|5.94|) and (|5.95|) . We set ^— i? = 



C and VL{i, R) = V{C, R) / R- Then ^Mli and (15:951) can be rewritten as 



V^^ + RVr = 0; i?>0, C + ^>0, (5.153) 
: + = 0; R>0, C + R = 0. (5.154) 

Changing variables from R to R = e^^ and letting t>{C,R) = V*{C,A), 



from (|5.153|) and (|5.154[) , for — oo < A < oo, we obtain 



Vl = Vl^ ■ C + e-^>0, (5.155) 
2e'^Vl+Vl+ [l + e-'^^)V* = Q ; ( + e'^ = 0. (5.156) 

We have thus reduced the problem to the heat equation, with A corre- 
sponding to the time variable, on a domain with a moving boundary = 
—e~^. On this boundary the condition (|5.156|) must be satisfied for all 
A G (—00,00). We could convert this problem to an integral equation, 
but this would not seem to have any advantages over (|3.48|) or (|5.148|) . 



5.5 Analysis near the corner {x,y) = (1,0) : 
x = l- 0(m-i/4)^ y = 0(m-i/2) 

Next we examine the behavior of 7r(fc, r) in a region that is further away 
from the boundary of ,t = 1, than the ^ scale. In this range the asymp- 
totics of -K{k, r) can be obtained by expanding the ray solution V{x, y) ~ 

81 



5/Cexp(^/5) as ^ (0,0). We set 



1 - X = 51/^ 



^, 1/ 



rV6 



and expand the ray solution V{x, ?/) as s — > and t ^ 00 letting 

t = -logs + T„ = To + 51/^X1 + 51/27^2 + 0(53/^), (5.157) 

s = 51/45,, 5. = 50 + 51/^51 + 51/252 + 0(53/^), (5.158) 

R = Y{S,,n;5) = Yo + 5^/^Y, + 6'^''Y2 + Oi6^/^), (5.159) 

r] = Z(5,, T,; 5) = Zo + 51/^^1 + 51/2^2 + 0(53/^). (5.160) 

On the {r], R) scale the ray expansion simplifies considerably, and with this 
simplified form we will be able to relate this scale to the (^, R) scale in 
subsection 5.4, as well as the (x, r) scale in subsection 5.2. In (|5.159|) and 
(|5.160|) it is understood that 5* and T* are replaced by their expansions in 
powers of 51/"^, as in (|5.157|) and (|5.158|) . Then we have F„ = F„(5o,To), 
Zn = Zn{So,To) and thus Yq = R, Zq = r] and F„ = Z„ = f or n > 0. 
Using (I5.157|l - (|5.160|) in (l330ll - (l335l) we obtain, after a lengthy calculation 



ii{k, r) 

where 

^(0) = 

^(1) = 



5ae-'^e^»/2 / ^ ^ /l + g-^o 
4e^° +2 + 77/ 



□—To 



-1/2 



exp 



^(0) ^(1) ■ 

1 h '^^ ' 

V6 51/^ 

(5.161) 



7/ / R 
2 V 1 + e-^o 

n2 



' 4 



i?log(l + e^o) 



1 



R 



2 V 1 + e-^o 



1 + e-'^o 2 V 1 + e-^o 



+ R 



R 



1 + e-^o 4 V 1 + e~^o 



R 



(5.162) 



3/2 



(5.163) 
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and 

-1 



R j '32' ' V 2 4 / V 1 + e-^o 



Here To is defined implicitly by 



+ | + ^e--ie-)]. (5.164) 



= 2To - ^ log5 - logi? + log (1 + e-^°) - 1. (5.165) 



Since Yq = R and Zq = rj imply that Sq = a/-R/(1 + e~^o) and = 
5o [2To - 2 log 5o - (log 5) /2 - 1] , we obtain (l5l65b . Note that Tq depends 



on and -R, and also weakly upon m, due to the term — log 5 = log m. 
To check that (|5.161|) satisfies (|5.5|) , we rewrite (|5.5|) in terms of rj and i? 



as 



(1 _ 51/4^ + 51/2^) (2 _ e*"/^^'') - (1 - e-*^/^'''-<5i/2i?e*-/^''' = 0. 

(5.166) 

Expanding \I' as 

^ = 51/2^(0) ^ ^3/4^(1) ^ ^^(2) ^ 0^^5/4^^ (5_^g7) 



and using this in (|5.166|) , we obtain at the first three orders (0(5^/^), 0(5^/^), 
0(5)) 



1 - e « 



(5.168) 



r/[v|/(o)]2_i?vl/(o) = 2vI/(0)vI/{i)+i?e*«'vl/«, (5.169) 
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and 



(1) 

R 



Differentiating both sides of (|5.165|) with respect to rj gives 

^1 



dTr 



V 



-To 



+ 1 + 



dv \2 y/R{l + e-To) 
From (|5.162|) we have 



1 + e-^o 



1 + e-^o 



R 



2 V 1 + e-^o 



+ 



4(l + e^of/2 2(l + e 



-To- 



+ 



R 



(5.171) 



5Tn 



2(l + e- 



drj ' 
(5.172) 



and thus, using (|5.171|) , we obtain 

^(0) = - 



R 



1 + e-^o ■ 

Similarly, differentiating both sides of (|5.165|l with respect to R gives 

-1 



(5.173) 



OR " 



-To 



1 + 



1 



From (|5.162|) and (|5.174[) we obtain 



^^-2V^ 



^(J) = - log (1 + e^o) 



(5.174) 



(5.175) 



We can easily check that (|5.173|) and (|5.175|) satisfy (|5.168|) . In the same 
manner, we can check that (|5.162|) - (|5.164|) satisfy (|5.169|) and (|5.170[) . 

We examine the asymptotic behavior of (|5.161|) as i? — with a fixed 
T] > (corresponding to So — > and Tq — > oo), which should match to 
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(|3.41|) when it is expanded as r ^ oo. Letting Tq — oo (and i? — 0) in 
(|5.165|) we obtain 

To ~ J log 5 + ^ log + ^ ("-^ + l] , (5.176) 



which when used in (|5.161|) gives 

, 5ae^'' , 1 f 7]^/R i? , , i?logi? r]^/R^ 



r]^VR rfR r]R^^^ 

~~8 ^ ^ 6 ^ T 



^ ^ m'-/2- exp (^-r/v^m^/^ - . (5.177) 

Expanding (|3.41|) as x ^ 1 and letting x = 1 — r]m~^/^ also lead to (|5.177[) , 
which verifies the matching between the {i], R) and (a;, r) scales. 

Next we expand (|5.161|) and (|5.165|) as ?7 ^ with a fixed R, by setting 
1] = ^5^^^^, which corresponds to To — oo and So with SqTq « 1. 
We note that in this limit, ^/R ^ oo. Then we obtain from (|5. 161)1 



2 



+ ri/V°/2i?V:R-:^(^2 + 5i/V^°/2X^ (5.178) 
and (|5. 165)1 becomes 

To ~ 5V4g-To/2^ + i log 5 + log i? + 1. (5.179) 
y R 2 

We let V = (5^/^e-^«/2^/v^ and using Tq ~ V + (log 6) /2 + log + 1 obtain 

4 ~ VeVe^/^ (5.180) 
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We also note that = \/5(^ / {RV'^) . Then we can rewrite (I5.178II as 

exp 



-ae I ( e e , r' 



i?V2v/2TV V 2V2 2V V 2(2 + V). 

(5.181) 



Since ^/R ^ oo, from (l5180b we obtain V ~ 2 log(^/i?) and (ISlSTl) thus 



agrees with (|5.122|) . This verifies the matching between the (r/, i?) and (^, i?) 



scales, for ,^ — >^ oo and 77 — > 0. 

5.6 Analysis near the corner (x, ?/) = (1,0): r = 0(1) 

We consider r = 0(1), which corresponds to there being only a few occu- 
pied secondary spaces, and x ~ 1. We shall discuss the scales k = m— 0(1), 



k = m — 0{^/rri), k = m — 0{^/mAogm) and k = m — 0{^/mlog'm) = 
m — y/m{\ogm)^. To get an idea of the forms of the expansions for n(k, r) 
on these scales, we first expand the results of subsection 5.5, which apply 
on the (r], R) scale, in the limit — 0, i? — 0, and then rewrite the results 
in terms of ^ and r. We note that 

l-x = b'lS = = -51/2 (log 5) e, R = S^^^r, (5.182) 

which relates the variables x, rj, ^ and ^. 



We rewrite (|5.165|) , which defines Tq, as 



^ ^ [2To-log5-logr + log(l + e-^°)-l] 



y/R \fr VI + e-^o 

(5.183) 



86 



Thus, if To = 0(1), 



Tn ~ - log ( |2 - 1 ) • 



(5.184) 



There is a singularity in this approximation to To when r = and this 
will lead to a transition in the asymptotics along r = We also note that 

To = if r = 2f , To < if r > 2^ , To > if f < r < 2^^ 



and (|5l83)l impl ies that if r < Tq oo with 



(5.185) 



We first consider the case r < Letting rj = m ^/^(logm)^, i? 
^-1/2^ and To — > oo in (|5.161|) yields 



.ae 



7r{k,r) ~ (5 — - — exp 



2 logm - rTo + - 



(5.186) 



As To ^ oo (151831) implies that 

log m log m log r 1 

To ~ = 1 \ — . 

° 2v^ 2 2 2- 



(5.187) 



which refines the approximation in (|5.185|) . Thus (|5.186|) becomes 

r logr 



7r(fc, r) 



~ exp 



log m + - log m 



ae 



-m 



r/2-fv^-l„-r/2 



(5.188) 



We note that the expression in (|5.188|) is algebraically small in m, and be- 
comes 0{m~^'/'^~^) when ^ = y/r. It is easy to see that (|5.188|) is precisely the 
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expansion of (|5.177[) as r] ^ {t] = m ^/^(logm) ^) and thus it also matches 
to the boundary layer solution (|3.41|) , when it is expanded as a; — >^ 1. 

When r > ^'^ we refine (|5.184|) by setting 



To = -log( ^-1 



logm 



Using (|5l89l) in BJM we obtain 

41og(r/e-0 + 2 



and thus 



To log I |2 - 1 



1-e/r 

^ 41og(r/e'-0 + 2 
(1 — ^'^/r) logm 



Using (151911) in (|5l611l leads to 

"2/2 C2 



n{k, r) ~ ae 



(5.189) 

(5.190) 
(5.191) 

(5.192) 



Next we examine the behavior of (|5. 161)1 in the transition region where 
^ ~ y^. We set 

To = log log m + log (5.193) 

and 



logm 



for 0, 4 ^ logm. Then we obtain from (|5.183|) 



4 ~ 2A/r log logm, 



and then defining ^* by 



log m / log m 



(5.194) 



(5.195) 



(5.196) 
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we find that (|5.183|) becomes 



^ = 21ogO - - logr - 1 + o(l). (5.197) 

This defines <) implicitly in terms of and r. Then on the (^*, r) scale we 
obtain from (|5.161)) 



n^k, r) ae " 



<> 



-m-''/2-l(log^)-2r^-r/2g-rv^_ 



(5.198) 



40 + 1 

As ^* ^ oo (I5l98l) should match to (l5J88b . From (I5l97b we see that 
^ oo and then using (|5.196|l , (|5.198|) becomes the same as (|5.188|) . We 
next check the matching between (|5.198|) and (|5.192|| . As ^ — oo, 0^0 
and (|5.196|) becomes 



P2 /i , 2e 41oglogm 
4^ ~ r IH h 



^yr log m log m 



(5.199) 



Using (I5l99|) in (l5192|l we obtain 

1/2 

exp 



7r(fc, r) 



ae 



2|* 



- log m ~ 2r log log m — log r 

2 2 



(5.200) 



Here we used —^* ^ log logm. The exponential part of (|5.200|) is the same 
as the exponential part of (|5.198|) . From (|5.197|| we obtain as ^ 



20' 



and thus 



<> 



rs^ \ / \ > r\j 







2|* 



1/2 



40 + 1 

Therefore, as ^* ^ -oo and ^ (I5l98l) matches to (I5l92l) . 



(5.201) 



(5.202) 



89 



Next we examine 7r{k,r) on the (^,r) scale by analyzing the balance 
equation (|2.3|) . We set 

1 



TT 



mylogm 



ir(^,r) exp[/(^) log m] 



(5.203) 



and rewrite (|2.3|) as 



2-_) ir(^,r)e^(«)'°s™ 
m/ 

m/ \ ymlogm 
H :r , — : 4 - 



r e 



/(?+l/(\A" log m)) log m 



/{C-l/(v^logm)) logm 



— ■ f g 

m — C,^/rn\ogm + r + 1 y ^/rnlogm^ 

+ J^/^ K{^, r + l)e^(«)i°g"^. (5.204) 

m — t,ym log m + r + 1 

From (|5.204|) to leading order (0(l/m)) we obtain the limiting equation 



r) [/'(O]' + r) - (r + r + 1) = 0. 



(5.205) 



By matching (|5.203|) to (|5.192|) , when r > ^'^ we must have f{^) 
Using this in (|5.205|) we obtain 

K{lr + l) r-e 



K{lr) 



r + 1 



Solving (|5.206|) gives 

"r-l 



no? 



and thus 



7r(A;, r) 



K{ll) T{r-e)K{ll) 
r! r(l-^2^ r! ' 



m A/log m r\ T{1 — ^'^) 
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-e/2. 



(5.206) 



(5.207) 



(5.208) 



We rewrite (I5.2U8I) as 

7r{k, r 



m^/\ogm r! ^"^^ r(l — 

We see that as ^ t V^> the factor T{r — ^'^) becomes singular, which indi- 
cates a transition in the asymptotics. Analysis of the range ^ ~ y/r will 
also determine H{^) in ^520% . and thus K{^,1) in BIM . To study the 
transition we scale ^ — = O [(logm)"^/^] with 



(5.209) 



i/logm 



, 7r(A;, r) = Hri®; m). 



Then (|2.3|) becomes 



+ 



2 7i:,(®;m) 

r + 1 ^ / log m 
h O ' 

m 



nr+i{®;m)+(l - ^) ® + 



m log m 
1 



m log m 



m 



(5.210) 



which we further approximate by 

(r + l)nr+i{®] m)-rnri®; m) + 



logm 



n''(®:m) + o 



logm 



0. (5.211) 



Before analyzing (|5.211[) we derive matching conditions for T-Cr{®',rn) as 
® ±oo. On the r) scale for f > y/r, (|5.188|) applies. This can also 
be obtained by analyzing (|2.3|) on the (^, r) scale for ^ > and using 



asymptotic matching. In either case we ultimately conclude that (|3.41|) 
remains valid for x ^ 1 as long as ,^ = (1 — x)y/m/ logm > ^/r. By setting 



C = a/^ + ®/-\/logm in (|5.188|) we obtain 



ae 



-m 



-r/2-l -r/2 -®V logm 



[1 + «(!)] 
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so that as ® +00 the matching condition for Hr{®; m) is 



+00. 



(5.212) 



By examining (|5.209|l as ^ — ^/r and using 



2v^ + 



v^logm 



-1 



vTogm 
(5.213) 



we obtain from (|5.209|) 



-2® rlv^r 



-00. (5.214) 



In view of the matching conditions we set 



Hri®; m) = m-'-/^^'e-®^^^[hr{®) + o(l)]. 



(5.215) 



With (ICTSb we see that 



(5.216) 



Also, since ^ ~ we have ^ < y/r + 1 and thus (|5.209[) may be used to 
approximate Hr+i{®; m) in (|5.211|) , as 



nr+ii®;m) 



m 



-r/2-1 



-®^r logm 



(r + 1)!' 



v^logm 

Using (|5.216|) and (|5.217|) in (|5.211[) we obtain to leading order 



(5.217) 



-®V2 



2v^n:(®) = ^H{V^) 



(5.218) 
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The solution to (|5.218|l that decays as 



-OO IS 



hr{®) 



g(v^) 



® 



(5.219) 



To determine H{y/r) we use the matching condition in (|5.212[) . Letting 
OO, we evaluate the integral in (|5.219|) as V2n and rewrite (|5.215|l , for 

OO, as 



r\\ r 



(5.220) 



Comparing (|5.212|) to (|5.220|) we obtain 



ae 



-r/2 



v^r(r + l^ 



so that 



ae 



(5.221) 



(5.222) 



Thus on the (®, r) transition scale we have 



7r{k,r) = Hr{®; m) 



ae~ 



m 



2V27r 



® 



-r/2-l^-r/2g-®v/Fk^ / C"" /^^M, r > 1, 

(5.223) 



and when ^ < y/r, from (|5.209|) and (|5.222|) , we obtain 

7r(A;, r) 



ae " m" 



r(i + f)r(r-a^-_,.^- 



r! 



(5.224) 



V27r mi/log m 

We note that letting ^, r — > oo in (|5.224|) and using Stirling's formula for 



the Gamma functions lead to (|5.192[) . This verifies the matching between 
the (^, r) scale and the ray expansion on the (x, y) scale. 
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We next examine r = 0. For ^ > 1 we can obtain the expansion of 
7r(/c, 0) as a limiting case of (|3.42|) , which leads to 



7r(fc,0)-(l-p)p^ 



ae 



-m 



(5.225) 



For < ^ < 1 we find that (|5.205|) still holds at r = 0, if we write 



7r(A:,0)-(l-p)p^ 



m 



my/logm 



Hence, 



ae m 



2/2 



271 m^/\og■m 



ei-«"r(i + ar(-e'), o<e<i. 



(5.226) 



We note that n{k, 0) — (1 — p)p'' is negative, since r(— < 0. There is a 
transition in the asymptotics when f = 1, and (|5.226|) also breaks down as 
^ — s> 0, due to the singularities of r(— ^^). We also note that on the ^ scale 



ae 



m 



a log m - ^ ( \og^ m 
1 + + ' ^ 



m 



(5.227) 



so that the geometric part of n(k, 0) is still asymptotically dominant for 
To study the transition range ,^ ~ 1 we let 7r(A;, 0) — {l — p)p^ ~ Ho(®; m) 



where now ® = ^/\ogm{^ ~ !)• Note that the transition point is the same 



for n{k, 0) and Tr{k, 1). By examining the matching condition to (|5.225|) as 
® +00 and to (|5.226|) as ® ^ — oo, we again set 



no{®;m) ~ m-3/2e-®Vi^;io(@) 



(5.228) 
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and from (15.21111 with r = we find that (logm)-^7Yo(®; m) ~ -Hi{®; m) 



so that /io(®) = — ^i(®) and hence 

7r(A;,0)-(l-p)p'= ~ _^^^-3/2g-®v/k^ /"° e-«V2^^^ @ = 

2\/2tx J-oo 

(5.229) 

We note that on the ^ scale the identity Xl^o — r, r) = (1 — p)p^ is 
asymptotically satisfied for ,^ < 1, in view of the fact that 

r=0 

Now we consider scales that have ^ — 0, so that m — k = o{^/m log m). 
If we define ^ by 



m 



k = \/ m log m ^, ^ = v^logm ^ 



then the factor m = exp(— logm/2) = e ^^/^ becomes 0(1) and 
(|5.209|) simplifies to 

vr(A;,r) ^ ^ e-^''/^ r > 1. (5.230) 

mlogm \/27i r 

Here we used T{r - f^) ^ (r - 1)! and H{0 ~ ae-^^/V^ir as ^ ^ 0. 
We have verified, after a lengthy calculation which we omit, that (|5.230[) 
can also be obtained by expanding m~^^'^^l{^, R) for ^ oo, R ^ with 
^ = 0{^/^\ogR), using the integral in (|3.47|) and our knowledge of 1^(0, R) 
(or i7(0, i?)) as i? ^ oo. If we expand 7r(/i;, 0) on the ^ scale for ^ ^ 0, we 
obtain from (|5.226|) 

^ 1 1 e-^'V2" 



7r(A;, 0) ~ ae'" 



m m 



'2^i 



(5.231) 
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so that on the ^ scale the geometric part of 7r(A;, 0) no longer dominates the 
right side of (|5.226|) . We thus re-examine the balance equation (|2.3|) along 



r = 0, i.e.. 



(2--]n{k,0) = (l--]n{k-l,0)+n{k + l,0) + -^n{k,l). (5.232) 
V m/ V m/ k + I 



On the ^ scale we let 7r(A;, 0) ~ m ^'v'(.^) and use (|5.230|) to approximate 
71 {k, I). Then (|5.232[) becomes asymptotically (after we multiply by m? log m) 



ae 



(5.233) 



The solution of (|5.233|) that matches, as ,^ ^ oo, to (|5.226)) is given by 

1 



9(0 = ae-^ 



(5.234) 



so that on the ^ scale with r = we have 

1 



7r(fc,0) 



ae 



(5.235) 



For ,^ ^ oo (|5.235|l reduces to (|5.231|| . This shows also that on the ^ scale 



{k = m — 0{\/m\ogm)) n{k, r) (r > 1) is smaller than 7r(/c, 0) by a factor of 
(logm)"^. 

Next we discuss the problem on the {n, r) scale where k = m — n and 



n = 0(1). Setting 7i{k, r) = C{n, r; m) we write (|2.3|) and (|2.4[) as 



a\ ^/ / Ci\ ^, m — n + \ 

2 C[n,r; m) = (1 C{n+l,r; m)H £(n— l,r; m) 

m/ V mJ m+r—n+1 



+ 



r + 1 



m + r — n + 1 



C{n,r + 1; m), n > 1 (5.236) 
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and 



(2 - —) £(0, r; m) = (1- —) £(1, r; m) + £(0, r + 1; m) 

\ m/ V m/ m + r + 1 



m + r + 1 

+ ( 1 - — ) £(0, r - 1; m), r > 1. (5.237) 



The comer condition (|2.6|) becomes 



2 - — ^ £(0, 0; m) = (^1 - £(1, 0; m) + ^^^^^(0, 1; m). (5.238) 



m + 1 



We expand C as 

£(n, r; m) = 



£(n, r) + £(^)(n, r) + 0(log" 

logm 



(5.239) 

The scale factor m~'^/^(logm)~^/^ must be included in view of matching 
considerations, which we discuss shortly. Using (|5.239|) in (|5.236|) - (|5.238[) 
we find that the leading term £(n, r) satisfies 

2L{n,r) - C{n + l,r) - C{n-l,r) = 0, r > 0, (5.240) 
2C{0, r) - C{1, r) - C{0, r - 1) = 0, r > 1, (5.241) 

with the corner condition 



2£(0,0) -£(1,0) =0. 



(5.242) 



The correction term C^^\n, r) also satisfies (|5.240|) - (|5.242|) . 

The general solution to (|5.240|) is C{n, r) = V(r) + nW(r) and we then 



obtain from (I5.241|l and (I5.242|l 



V(r) - V(r - 1) = >V(r), r > 1, 



(5.243) 
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and 



V(0) = W(0), (5.244) 

which implies that C{n, 0) = (ra + l)V(O). A similar argument can be used 
to determine £^^)(n, r) as V^^\r) + nW*^^''(r) so we write the expansion on 
the in, r) scale as 



C(n,r; m) 



m^/^i/logm 
1 



+ 



logm 



{ V(0) + ^W(/) + r2W(r) 
1=1 

r 

v« (0) + 5^ >v« (0 + ^w^'H^ 



} (5.245) 



for r > 1, and 

£(n, 0; m) 



n + l 



3/2^1^ 



V(0) + 



logm 



(5.246) 



Now we try to match (15.24511 and (I5.246D to (I5.230D and (15.2351) , noting 



that n = iy/m logm. On the ^ scale £(n,0) in (|5.246|) becomes O 



m 



which is of the same order as (|5.235|) , but the expansions cannot match 
since (|5.246|) will be linear in n (or while (|5.235|) does not vanish as ,^ — *• 
0. Problems also arise in matching the n and ^ scales for r > 1. For a fixed r 
and n oo, the leading term in (|5.245|) becomes m"'^/^(logm)^^/^n>V(r) = 
m~^^>V(r) which is larger than (|5.230|) by a factor of logm. We must thus 



set >V(r) = 0. We can match (|5.230|) to the correction term in (|5.245|) by 
setting W(^)(r) = ae-" / {^/2^lr). Then we would have, for r > 1, 



C{n,r; m) 



m' 



3/2 



m 



V(0) + 



logm 



n 
r 



E 



ae 



+ v(^)(o) 



(5.247) 
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and 

T). -1-1 r 1 ' 

(5.248) 



L[n, U; m) 



v(o) + -^v(^)(o) 

logm 



We have verified, by numerical computations, that the order of magni- 
tude of 7r(/c,r) on the {n,r) scale does seemtobe 0(m~^/^(logm)~^/^), that 
7r(/c, 0) is approximately proportional to + 1, and that Tr{k, r) is approxi- 
mately constant for r > 1 and k = m — 0{1). However, the problems with 
the matching suggest that there is yet another scale in the problem, which 
corresponds to n — > oo and ,^ ^ 0. We have not been able to identify this 
new scale. 

To summarize this subsection, we obtained results for n{k,r) on the 
r) scale, treating separately the cases ^ > y/r, ^ ~ y/r, and ^ < y/r, 
for r > 1. For r = we gave results for ,^ > 1, ,^ ~ 1, and ,^ < 1. For 
^ = ^ -\/log m = 0(1) we obtained the simplified result in (|5.230|) for n(k,r) 
for r > 1. For 7r(fc, 0), (|5.235|) applies on the ^ scale. On the {n, r) scale we 
obtained (|5.247|) for r > 1 and (|5.248|) for r = 0. However, there is still a 



"gap" in the asymptotics between the n and ^ scale, a gap which we have 
not been able to fill. 

6 Numerical Studies 

We assess the accuracy of some of the asymptotic formulas we obtained. 

In Table 1 and Table 2 we test our asymptotic results for m ^ oo with 
a fixed p < 1. Table 1 has 7r(0, 1), where (|3.16|) applies with k = and 
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r = 1, 7r(m/2, 1), where (|3.14|) applies with x = 1/2 and r = 1, and 7r(m, 0), 
where (|3.20|) applies with n = 0. We consider p = 0.5 and increase m from 
10 to 30. The agreement is not particularly good for {k, r) = (0, 1) but does 
improve significantly as m increases. For {k, r) = {m/2, 1) (boundary layer 
near ?/ = in subsection 4.2) and (k, r) = (m, 0) (comer layer in subsection 
4.3) the agreement is good even for m = 10, with errors of at most 10%. 
In Table 4, we consider {k,r) = (0,m), (m/2,m), {m,m), so that y = 1 
and X =0, 1/2, 1. The asymptotic formulas that apply are now (|3.9|) (with 
(I3l0l) - (l3l3l) ), (ED (with - (I321»), and (with x = y = 1). We again 
consider p = 0.5 and m = 10, 20, 30. Now we obtain generally excellent 
agreement, with the worst error in Table 2 being about 3%, which occurs 
for A; = and r = m = 10. 

These comparisons show that the asymptotics agree reasonably well 
with the exact numerical values of n(k, r). The comparisons also clearly 
demonstrate the necessity of analyzing separately the different ranges of 
(/c, r), when m oo. 



100 



Table 1 

0.5; (A;,r) = (0,1), (m/2,1), (m, 0) ; 10 < m < 30. 



{k,r) 


m 


exact 


asymptotic 


(0,1) 


10 


5.83 X 10-6 


9.76 X 10-6 




20 


1.78 X 10-*^ 


2.38 X lO-'^ 




30 


8.43 X 10-13 


1.03 X 10-12 


(m/2,1) 


10 


8.57 X 10-5 


7.81 X 10-5 




20 


7.92 X 10"* 


7.62 X 10-^ 




30 


7.61 X 10-11 


7.45 X 10-11 


(m, 0) 


10 


2.71 X 10-^ 


2.44 X 10-^ 




20 


2.51 X 10-^ 


2.38 X 10-^ 




30 


2.11 X IQ-i" 


2.32 X IQ-io 



Table 2 



0.5; (/c, r) = (0, m), (m/2,m), (m, m); 10 < m < 30. 



(/c,r) 


m 


exact 


asymptotic 


(0,m) 


10 


7.59 X 10-13 


7.40 X 10-13 




20 


3.61 X 10-2' 


3.59 X 10-25 




30 


1.52 X 10--^^ 


1.51 X 10-37 


(m/2, m) 


10 


4.86 X 10-9 


4.81 X 10-9 




20 


4.93 X 10-1^ 


4.90 X 10-17 




30 


5.02 X 10-25 


5.00 X 10-25 


(m, m) 


10 


3.67 X 10-^ 


3.57 X 10-7 




20 


3.45 X 10-13 


3.41 X 10-13 




30 


3.28 X 10-19 


3.25 X 10-19 



101 



References 



[1] Aldous, D., Some interesting processes arising as heavy traffic lim- 
its in an M/M/oo storage process. Stochastic Processes and their 
Applications 22 (1986) 291-313. 

[2] Coffman, E. G., Jr., Flatto, L. and Leighton, F. T., First fit allocation 
of queues: Tight probabilistic bounds on wasted space. Stochastic 
Processes and their Applications 36 (1990) 311-330. 

[3] Coffman, E. G., Jr., Kadota, T. T. and Shepp, L. A., A stochastic model 
of fragmentation in dynamic storage allocation, SI AM J. Comput. 14 
(1985) 416-425. 

[4] Coffman, E. G., Jr. and Leighton, F. T., A provably efficient algorithm 
for dynamic storage allocation, J. Comput. System Sci. 8 (1989) 2-35. 

[5] Coffman, E. G., Jr. and Mitrani, I., Storage of the single-server queue, in 
Queueing Theory and its Applications, Liber Amicorum for J. W. Cohen, 
edited by O. J. Boxma and R. Syski, CWI Monogr., 7, North-Holland, 
Amsterdam (1988) 193-205. 

[6] Courant, R., and Hilbert, D., Methods of Mathematical Physics, Vol. 2, 
Wiley - VCH, Wiley Classics Libray 20 (1989). 

[7] Knessl, C, Asymptotic expansions for a stochastic model of queue 

storage. Annals of Applied Probability 10 (2000) 592-615. 

102 



[8] Knessl, C, Geometrical optics and models of computer memory frag- 
mentation. Studies in Applied Mathematics 111 (2003) 185-238. 

[9] Knessl, C, Some asymptotic results for the M/M/oo queue with 
ranked servers, Queueing Systems 47 (2004) 201-250. 

[10] Kosten, L., Uber Sperrungswahrscheinlichkeiten bei Staffelschaltun- 
gen, Electra Nachrichten — Technik 14 (1937) 5-12. 

[11] Newell, G. E, The MjMj oo Service System with Ranked Servers in Heavy 
Traffic, Springer, New York (1984). 

[12] Preater, J., A perpetuity and the M/M/oo ranked server system. 
Journal of Applied Probability 34 (1997) 508-513. 

[13] Sohn, E. and Knessl, C., A simple direct solution to a storage alloca- 
tion model. Applied Mathematics Letters 21 (2008) 172-175. 

[14] Sohn, E. and Knessl, C., The distribution of wasted spaces in the 
M/M/oo queue with ranked servers. Advances in Applied Probability 
40 (2008) 835-855. 

[15] Sohn, E. and Knessl, C., Storage allocation under processor sharing I: 
Exact solutions and asymptotics, submitted for publication. 



103 



